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1.0  INTRODUCTION 


This  document  contains  a  bibliography  o£  reports,  papers  and  books 
that  deal  directly  with  Prony's  method.  No  attempt  has  ben  made  to 
Include  papers  which  fall  under  the  general  categories  of  approximation 
with  exponential  sums*  linear  prediction,  parameter  estimation  or  system 
Identification  even  though  these  topics  are  closely  related  to  Prony's 
method.  The  papers  which  we  found  useful  In  these  fields  are  referenced 
In  the  other  volumes  which  make  up  this  report. 

Section  2.0  of  this  volume  contains  a  translation  from  the  French  of 
Prony's  original  paper.  Ihls  has  been  Included  for  two  reasons.  The 
first  being  that  while  Prony's  original  paper  Is  highly  referenced,  we 
suspect  that  most  people  have  never  seen  It  let  alone  have  read  It.  The 
second  reason  for  including  it  Is  that  we  find  It  extremely  Interesting 
to  read. 

Section  3.0  contains  an  annotated  bibliography  of  Journal  articles, 
reports  and  theses  that  deal  with  Prony's  method.  In  most  cases,  the 
annotation  was  taken  from  the  author's  abstract  of  his  paper.  Section 
4.0  contains  lists  of  papers  presented  at  conferences,  symposia  and 
meetings.  Section  5.0  contains  a  list  of  books  dealing  with  the  subject 
of  Prony's  method. 

It  Is  hoped  that  the  usei's  of  this  list  will  cross  reference  It 
with  their  own  bibliographies  on  Prony's  method  and  will  let  ETI  know 
of  any  additions  which  they  find. 
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2.0  TRANSLATION  OF  PEONY’S  ORIGINAL  PAPER 


EXPERIMENTAL  AND  ANALYTICAL  ESSAY 
ON 

THE  EXPANSION  PROPERTIES  OF  ELASTIC  FLUIDS 
AND  ON  THE  FORCE  OF  EXPANSION  OF 
WATER  VAPOR  AND  ALCOHOL  VAPOR  AT  DIFFERENT 
TEMPERATURES 

By  R.  Prony,  1795 
Tranilated  by  Dr.  Ann  Sanders /ETI 

GENERAL  CONSIDERATIONS 

In  approximately  the  last  40  years »  Physics  has  been  enriched  greatly 
by  a  large  number  of  observations  made  with  extrema  care  by  clever  and  ex¬ 
perienced  men.  This  knowledge  Is  Increasing  dally,  and  Its  substancs  Is 
becoming  ever  more  precious  as  the  perfection  of  new  Instruments  provides 
more  and  more  precision  to  the  experiments;  phenomena  have  bean  correlated, 
compared  and  classified;  the  language  of  an  Important  part  of  science  has 
become  analytical;  reasonable  theories  have  replaced  futile  and  often 
absurd  theories,  which  had  bean  caught  until  the  middle  of  this  century. 

The  study  of  nature  In  terms  of  Its  observation  and  the  knowledge  of 
Its  operation,  seems  to  me  to  consist  of  two  aspects  which  should  not  be 
confused:  the  explanation  of  effects  and  their  measurement . 

The  explanation  of  effects  consists  of  findings.  In  a  class  of  multi¬ 
faceted  phenomena,  those  simple  or  basic  phenomena  from  which  the  others 
are  derived,  or  of  which  they  consist  In  various  combinations.  It  also 
consists  of  showing  how,  from  the  most  varied  appearance,  one  can  untangle 
the  nporatlon  and  the  being  of  the  elements  considered  to  be  the  basis 
of  the  system. 
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Thus,  by  starting  from  the  affinity  of  certain  substance?  which  are 
conulv'iered  as  basic  phenomena.  It  has  been  found  that  meteorological 
phenomena,  combustion  phenomena  etc.,  result  from  these  same  affinities, 
and  show  up  under  different  disguises  which,  until  recently,  hid  their 
nature  from  physicists.  It  Is  by  this  decomposition  of  complex  effects 
Into  basic  effects  that  we  can  discover  certain  secrets  of  nature,  which, 
while  It  allows  us  to  lift  one  comer  of  the  veil  which  covers  It,  keeps 
another  comer  tied  by  a  knot  which  our  hand  cannot  untie. 

The  measurement  of  effects  is  the  evaluation  of  the  degree  to  which 
each  effect  Is  varlabJe  when  one  varies  either  the  causes  which  produce 
It,  or  other  effects  to  which  It  is  tied.  It  Is  known,  for  Instance,' 
that  the  tendency  of  fluldo  to  vaporise  Is  on  the  one  hand  enhanced  by 
temperature,  and  on  the  other  hand  reduced  by  the  atmospheric  pressure, 
and  that  vaporization  occurs  only  when  the  former  force  exceeds  the 
latter  one;  In  addition  to  this  fact.  It  may  be  necessary  to  give  the 
values  of  the  pressures  which  produce  vaporization  equilibrium  at  different 
temperatures.  This  same  kind  of  reasoning  is  applicable  to  Infinitely 
many  other  phenomena. 

It  can  thus  be  seen  that  the  explanation  of  effect'!,  whose  great 
advantage  is  to  simplify  science  and  to  correlate  its  various  aspects, 
by  analysis  and  by  decomposition  Into  basic  phenomena,  must  be  complemented 
by  the  measurement  of  these  same  affects,  which  is  always  very  useful  and 
often  Indispensable  when  It  is  desired  to  apply  theoretical  discoveries 
to  the  needs  of  society. 

Only  experiments  can  provide  the  Initial  data  on  the  measurement  of 
physical  effects;  following  those,  however,  calculations  can  he  very  use¬ 
ful,  either  to  obtain  results  Intermediate  to  those  obtained  from  the 
experiments,  or  to  correct  experimental  anomalies.  The  method  used  in 
that  case  is  called  interpolation.  Its  aim  Is  to  find  a  relatlonohlp 
between  two  or  three  variables,  such  that  when  a  certain  value  is  given 
to  one  or  two  of  the  variables,  a  equally  specific  value  la  obtained  for 
the  second  or  third  variable  when  the  problem  is  considered  In  that  manner. 
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it  can  be  solved  In  an  infinite  number  of  ways,  since  there  are  an  Infinite 
number  of  functions  which  can  be  obtained  by  the  same  substitutions.  It 
would  be  erroneous,  however,  to  think  that  all  those  solutions  are  equally 
applicable  to  a  given  case.  Even  though  nature  is  governed  by  seemingly 
few  and  simple  general  laws.  It  has  as  many  modifications  in  its  workings 
as  variety  in  its  forms,  and  each  measurable  phenomenon  is  always  corre¬ 
lated  to  a  specific  function  which  represents  It  exclusively. 

The  problem  of  Interpolation  thus  has  two  distinct  parts;  the  first 
part  consists  of  satisfying  given  numbers,  and  the  second  part  consists 
of  finding  the  function  applicable  to  the  phenomenon  amongst  all  the 
functions  which  satisfy  the  conditions. 

In  ^19  of  my  Analysis  Lessons,  I  have  given  a  solution  to  the  first 
part  of  the  problem  which  Is  used  very  often,  especially  as  a  method 
of  correction.  On  the  same  subject;  Lagrange  has  published  a  beautiful 
treatise*  where  he  handled  the  subject  in  a  more  general  way  than  had 
been  done  before.  The  students  which  know  the  theory  explained  in  my 
lessons  if  19,  20  and  21,  will  be  able  to  study  this  work  with  ease, 
and  will  be  able  to  gain  much  from  the  time  they  will  have  spent  on  it. 

From  the  knowledge  we  have  now,  the  solution  to  the  second  part 
does  not  seem  to  be  governed  by  general  rules,  especially  whan  the  number 
of  observations  is  small,  and  when  the  observations  do  not  span  a 
wide  range.  The  only  guidelines  available  in  this  painful  search  are 
careful  attention  to  all  the  details  and  the  progression  of  the  experl- 
mentS)  repetition  of  measurements,  and  analogies.  These  difficulties, 
combined  with  the  difficulty  of  precision  in  experiments,  are  the 
reason  why  exact  determinations  of  physical  laws  are  very  rare  in  Physics. 

In  1790,  I  had  the  opportunity  to  follow  very  detailed  and  careful 
experiments  on  the  expansion  force  of  water  vapor,  and  it  was  my  aim 
to  find  the  applicable  formula.  The  rngularlty  of  the  data  led  me  to 
believe  that  the  task  would  be  easier  than  it  actually  was.  However, 
after  some  effort,  I  found  a  certain  function  which  not  only  perfactly 
expressed  the  relationship  between  the  temperature  and  the  spring  of  the 


*  Footnote  Illegible. 
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watery  vapor,  but  which  seemed  to  apply  in  general  to  phenomena  depending 
on  elastic  fluids.  I  applied  them  to  careful  experiments  which  Prleur 
made  on  the  expansion  properties  of  air  and  of  different  air-form  fluids; 
this  strengthened  my  opinion,  and  I  decided  to  publish  my  results. 

The  first  observation  which  directed  me  to  the  true  form  of  the 
function  was  the  consideration  of  certain  geometric  progressions  which 
are  present  In  certain  phenomena  related  to  elastic  fluids.  One  of  the 
most  remarkable  examples  of  this  Is  the  relationship  between  the  density 
of  atmospheric  layers  and  their  respective  height.  Since  this  relation¬ 
ship  Is  exponertlal,  I  suspected  that  In  other  cases,  where  one  such 
quantity  would  not  be  sufficient,  two  or  more  could  be  Introduced.  This 

XXX  X 

led  to  an  equation  of  the  form  t  ■  w,p,  +  +  M,p,  +  ...  +  y  p 

i  i  nn 

where  t  and  n  are  the  two  variables  and  *^2“’ 

constants  appropriate  to  the  given  phenomenon. 

It  Is  known  that  the  above  equation  results  from  an  equation  with 
linear  finite  differences  or  gives  the  general  term  in  a  recurrent  series 
of  order  n.  However,  series  of  the  type  where  a  given  term  la  obtainable 
from  a  certain  number  of  preceding  terms  seems  Indeed  to  be  compatible 
with  natural  effects  where  elasticity  plays  a  large  part.  The  conserva¬ 
tion  of  living  forces  which  is  contained  In  this  property  of  bodies  always 
causes  the  actual  state  to  be  dependent  on  the  previous  states.  The  re¬ 
search  of  Lagrange,  which  I  mentioned  above.  Is  also  based  on  recurrent 
series.  He  has  given  several  methods  to  find  those  which  must  interpolate 
a  given  series,  where  one  can  observe  the  elegance  and  the  depth  which 
one  comes  to  expect  from  such  a  great  analyst.  Since  the  method  I  have 
used  In  my  calculations  differs  from  his  methods,  which  I  did  not  know 
when  I  started  this  work,  I  will  now  describe  It. 
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Interpolation  Method  Apollcable  to  Phenomena  PependiEg  on  Elastic  Fluids 


The  experiments  must  be  directed  as  much  as  possible  cowards  yielding 
equidistant  results;  when  this  condition  cannot  be  met  (which  must  happen 
very  rarely) »  but  when  the  results  are  nevertheless  numerous  and  close 
enough  together,  they  can  be  made  to  be  equidistant.  This  can  be  done 
either  by  graphic  methods,  by  drawing  the  experimental  curve,  or  by  cal¬ 
culations,  namely  by  considering  three  consecutive  results  Z  ,  Z  and 

1  ^11 

^111*  second  and  the  third  are  a  distance  x  and  x  respec¬ 

tively  from  the  first.  The  result  Z,  which  belongs  to  Che  series  of 

equidistant  points,  and  which  Is  a  distance  x  from  Z. ,  can  be  calculated 

^  th 

from  the  following  equation,  which  Is  obtained  from  those  In  my  19 
lesson  of  analysis: 


‘ill 


) 


This  formula  can  be  simplified  greatly  by  calculating  the  difference 
between  Z  and  Z^.  By  setting  Z^^  -  Z^  ■  and  "  ^1  " 
obtains 


z  -  Zl  “  — 


X  -  X 


,x  -  X  1  X  -  X  iiv 
( - 1 —  U  -  U  )  . 


In  order  to  Ld  mistakes.  It  Is  advisable  to  apply  both  the  formula 
and  the  graph::.c  •aethod,  which,  when  applied  carefully  and  on  a  large 
scale,  will  give  a  precision  comparable  to  Che  experiments  themselves. 


Once  this  preparation  has  been  made  (and  as  I  said,  the  cases  where 
this  Is  necessary  will  be  extremely  rare) ,  a  certain  number  of  equidistant 
results  will  be  taken  which  will  encompass  either  the  entire  experimental 
range,  or  a  large  part  thereof.  Then  the  variable  Z  Is  taken  which 
indicates  the  measurement  of  successive  effects  which  correspond  Co  cer¬ 
tain  values  of  another  variable  x,  which  Indicates  to  which  term  on  a 
given  scale  the  effects  Z  approach. 


To  satisfy  a  number  2n  of  results,  one  will  generally  have 


•  *  “(n)‘’(n) 


and  to  satisfy  a  number  2n  +  1  of  results, 

^  “  “1^  *  •  •  *  "(nAn)  *  +  1) . 

Where  u^,  .  etc.,  and  p^,  ®tc.,  are  constants  whose 

values  are  determined  from  the  experdmental  results,  as  will  soon  be  seen. 


I  have  given  two  general  formulas,  although  usually  one  of  the  two 
would  be  sufficient,  but  I  had  a  simplification  In  mind  which  is  an 
important  Inclusion  in  the  method,  and  which  will  now  be  described.  The 
determination  of  p^,  p^^,  etc.,  depends  on  the  solution  of  an 

equation,  and  by  using  the  second  formula,  one  satisfies  an  odd  number 
of  observations  by  an  equation  which  Is  of  the  same  order  as  that  required 
for  satisfying  the  next  smaller  even  number  of  observations.  Thus,  one 
can  satisfy  four  or  five  observations  by  calculating  a  second  order 
equation,  six  or  seven  observations  by  calculating  a  third  order  equation, 
eight  or  nine  observations  by  calculating  a  fourth  order  equation,  and 
so  on.  It  will  almost  never  happen  that  eight  or  nine  results  oust  be 
included  in  the  formula,  and  without  exceeding  the  limits  where 
methods  for  the  solution  of  numerical  equations  exist,  one  will  be  able 
to  treat  all  the  cases  which  Physics  ordinarily  presents.  In  addition 
to  that  advantage,  it  should  be  mentioned  that  there  will  be  at  most  only 
half  as  many  values  of  Z  as  there  are  observation  points.  Instead  of  with 
the  formulas  which  refer  to  parabolic  curves,  and  which  always  have  as 
many  terms  as  there  are  observation  points. 


Here  then  Is  the  method  by  which  the  constants  in  equations  (1) 
and  (2)  can  be  determined,  using  the  given  results. 


First  case,  whera  the  nuabar  of  obsarvatlons  Is  avan: 


In  my  analysis  lassons  numbars  20  and  21,  I  hava  shown  Chat  aquation 
(1)  producad  tha  ganaral  Cam  In  a  racurrlng  sarlas  of  ordar  n,  and  also 
that  tha  tarns  In  such  a  sarlas,  takan  at  cartaln  aqual  Intarvala,  always 
producad  sarlas  of  tha  sama  ordar.  Glvan  that,  lat  thara  ba  tha  following 
two  sarlas,  tha  first  of  which  producaa  tha  obsarvad  rasults  of  tha 
particular  valuaa  of  Z  fumlshad  by  tha  axparlmant  and  tha  aacond  of  which 
consists  of  tha  corraspondlng  valuaa  of  x,  namaly  tha  obaarvad  rasults 


Zq;  Z^j  2(^2) ••••^Zn-l 


and  tha  corraspondlng  valuas  of  x 


0;  x^;  2x^,...nx^;  (n+l)x^;  (n+2)x^....(2n-l)x^. 


Tha  quantltlaa  Z^,  Z^,  Z^^,  ate.  must  form  a  racurrlng  sarlas  whosa  ordar 
must  ba  found. 


Lat  A^t  A^,  A^^  ba  undataralnod  coafflclants,  satisfying  tha  ralatlon- 

ships 


Zq  +  Zj  +  Aj,  Z^^  + . -f  Z(„,  -  0 

h  *  \  hi  *  * . *  '^rn) 

^11  *  \  ^111  *11  ^Iv  . *(n) 

'ill  *  *1  'iv  *  *11  ♦ . *  '(n)  '(nt3) 


0 

0 

0 


An  Z,  "^A.  Z ,  .  ■^■A  Z,  + 
0  (n-1)  1  (n)  11  (n+1) 


*(n) 


(For  convanlence,  A^^^  may  ba  sat  to  1  In  tha  numarlcal  applications.) 


These  a  equations  will  produce  u  ratios  -v—  '  ‘  '  *  ’  *  ’  4  ~ 

^*(n)  ^(n)  ^(n)  ^(n) 

which  mahe  up  the  desired  order  of  the  series,  and  one  will  obtain: 


For  n  ■  1 

A(,  :  Ai  -  -2i  :  Zg 

Value  obtained  from  two  observations. 


For  n  “  2 


*0  •  u  z,z^^  -  z^  z^ 

A  .  A  -  '  V iil 

^  "0^11  - 

Values  obtained  from  four  observations. 


For  n  ■  3 


A-  •  A 


.  <^Ai  - 

111  -  ZiiZiv>Zo  *  ■ 


t  tv  11  Iti'  1 


11  11 


iAliifv 

^iW^ll 


*1  ‘  *111 


*<ZlZlv  -  ^Ii^iii’^iu  *  <^11^11  -  Vlv^^lv  *  *^0^111  ~  ’^l^ll^  \ 


-^hhy  -  hi^Ui>h  -  <"ll"ll  -  Vlv>^l  -  <"o"lll  -  "l"ll)  ^iv 


.  ,  -  ^1^111>^1:  .  *  <"0^111  -  ^^l>"lv  ^  «l"l  -  Vll»v 

*11  =  *111  ■  -(Z^^Z^^  -  Z^Z^^^)Zj^-  -  (Z„Z^j^  -  (Z^Z, 


Values  ('  stained  from  six  observations. 


Then,  upon  solving  the  equation. 


+  A.  a  +  A. .  a  +  A.,,o  + . ...+A,.o  -  0, 
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i  1  i  1 

the  n  roots  which  will  be  obtained  will  be  the  n  values  p.  ,  p._,  P..j***P/  \ 

1  11  111  vn; 

from  which  the  values  of  P^,  P^^,  ®tc.  will  be  obtained.  Finally,  the 

quantities  «tc.  will  be  given  by  the  equations: 


.  «  -  . »  - 

1  X  X  X.  X  X  X  X,  X, 

<»i  -"u'  '“l  . -‘“i 

(z  -  (z  -  . (z  - 

^ii  X  X  X  X  X  X  X.  X. 

^‘^ii  "  “^i  ^  ^^11  “  ^iii^  ^^il  "  ‘’iv^ . ^^il  “  ^(n)^ 

(Z  ■  o})  (Z  ^  p^{)  (Z  -  p"^) . (Z  - 

^ili  X,  X.  X  X  x^  X  X  X 

^‘^ili  “  ^i  '  ^‘^iii  "  ‘’ii^  ^‘’iii  “  “^iv^ . ^^ill  “  “^(n)^ 


t 


I 


I 


9 


‘(n) 


X  X  X 

(Z  -  p^^)  (Z  -  p^^)  (Z  « 


(p 


(n) 


-  "i'*  ‘"(i)  - 


‘’ii>  ‘"(i)  ■ 


"iii^ 


(p 


i  - 

(n) 


by  noting  that  in  the  numerators  all  the  coefficients  of  the  powers  of  Z 

must  be  changed  to  subscripts  of  the  same  order,  i.e.,  that  Z^  must  be 

replaced  by  Z  (or  that  all  the  terms  which  do  not  contain  Z  must  be 
0  2 
multiplied  by  Z^) ,  that  Z  must  be  replaced  by  Z^ ,  .hat  Z  must  be  replaced 

by  Z^^,  etc.,  thus  one  will  have  in  the  case  of 

n  -  l...w^  -  Zq 

To  satisfy  two  observations. 
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To  iiatlafy  six  obsorvatlons 


This  agrees  with  the  formulas  which  I  gave  In  my  20th  lesson  of  analysis 
when  n  Is  set  to  1.  In  order  to  make  t*  elimination  calculations  which 
give  u  ,  u..t  etc..  Identical  to  those  p.usented  In  the  lesson,  one  can 
*1  *1 

set  "  'I'i  ’  *^11  "  o"®  can  look  for  etc..  In 

terms  of  i|»^,  etc.  One  knows  that  here  n  Is  thy*)  constant  Increase 

in  X,  or  Ax.  Once  the  quantities  p^,  «tc.  and  ®l*c. 

are  determined,  their  values  can  be  substituted  the  equation  (1). 


'(n)  "(n) 


which  will  then  be  dlapoaed  to  aetiafy  the  2n  given  obaervationa*  and 
which  will  be  able  to  produce  aome  reault  Intermediate  between  thoae 
obtained  in  fact. 

Second  caae,  where  the  number  of  obaervationa  ia  odd: 

In  order  to  aolve  the  aecond  caae*  it  ahould  be  noted  that  aquation 
(2)  differa  only  by  the  conatant  aquation  (1).  Thua  the  aariea 

of  u  obtained  from  (2)  ia  of  the  aame  nature  aa  that  obtained  from  (1),  the 
only  difference  being  that  in  (2)  each  term  ia  increaaad  by  If 

thaaa  terma  are  therefore  decreaaed  by  u  >  ...»  the  remaindara  will  aatiaf y 
the  relationahipa  which  make  up  a  aariea  of  order  n*.  Thua»  by  preaerving 
the  notation  of  the  previoua  article »  one  haa 

^0^V^(n+l)^  ”*■  ^ii^^ii'^^Ctri-l)^  ■*‘*-*'*‘  ^(n)^^(n)“‘^(n+l)^  *  ° 

^0^^i“^(n+l)^  *'■  ^i^^ii“’^(n*fl)^  *'■  ^ii^^lii"’^(n+l)  ^(n)  ^^(nfl)"‘'(nfl)^  “  ° 

^0^^1^^Itfl)^  ^l^^lii“^n4-l)^  *  ^ii^^iv’^n+1)^  ^n)^^(IVf2)"^rrt•l)^  "  ° 

^O^^iii'^'^n+l)^  \^^lv“‘^(n+l)^  *’■  ^ii^V’^Cnfl)^  ^(n) ^^tri-3)'*’^(nfl)^  "  ° 


hi 


A„(Z 


0'“(n-l)"‘^(n+l)^  ^i^^(n)“‘^(n+l)^  ^ii^^(rH-l)“*^(nfl)^  •"'‘^(n)  ^^(2n-l)“'^(n+l)^“°  ^ 


)  +  Aj 

,1 

'0 

If  theae  equatioua  are  subtracted  from  one  another*  u 


^0^^(n)“‘*(n+l)^  \^^(n+l)^"^(n+l)^  ^il^^n+2)^“^(n+l)^  ■*■*  * ^^(2n)“‘^(tH-l)^“° 


eliminated,  and  the  equations  will  become 


(n+1) 


will  bo 


Vi 


*  The  recurring  aeries  ia,  in  fact,  of  order  n-t-l,  but  the  equation  of  the 
relationship  haa  a  root  equal  to  unity:  i.e.,  in  the  term  ^(nfl) 

one  has  p  . .  ■  1,  such  that  the  order  of  the  equation  may  be  decreased 
by  one. 
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*0  ^  *0  *  *1  “l  *11  *  "ll  *••••*  *(n)  ‘  "(n)  •  ° 

*0  *  ^1  *  *1  *  ’^11  *  *11  *  ^111  *••••*  *(n)  *  ^n+l)  ■  ° 
^  *  ^1  *  *1  *  ^111  *  *11  *  'iv  *••••♦  *(«)  *  ^n+2)  ■  “ 
*0  *  'ill  *  *1  *  'iv  *  *11  ^  ^  \„)  ^«.3)  -0 


*0  *  '(n-1)  *  *1  *  '(n)  *  *11  *  '(iH-1)  *•■•■*  *U)  *  '(2"-l)  ■  “ 

For  tho  ookt  ot  convonlonca,  ono  con  ••c  A  ^  ■  1  In  tha  nuaarloal 
application! a 

^0  ^11 

From  tho!!  aquatlona,  the  valuas  *  a  *  •tc#,  can  ba 

^(n)  ^(n)  ^(n> 

obtalnad,  namaly: 


For  n  •  la  a  a  Aq  t  ■  A  i  A  Zq 


Valua  obtalnad  from  thraa  obiarvatlonia 


For  a  ■  2a  a 


A  Z.  L  Z,,,  -  til,,  L  I, 


A  a  A  ■  *  ***  **  ■  -*■* 

*0  •  *11  1  Zg  *  'll  '  *  'l  *  'l 


*1  ‘  *11 


A  Z,  A  -  1  Zg  A  2,,, 

*  'o  *  'll  -  ‘  'l  *  'l 

Valuai  obtained  from  flva  obiarvatlonia 


For  n  ■  3a  a ■ 


_  ^(l^llZZ,,,-AZ,,A2,^)AZ,  a-(Ai^AZ,^-A2,,AA  )AA,„MA2,,A2^,-AZ^AZ  )A2, 

0  til  -AZiiiA2;^;-Az“AZ;,)AZg-?AZ^Ai^^;AZjiM^j;5AZ^-(AZ^jAZjj-AZ^AZ^^^)AZjj 
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^4  *^«lll  m  • 

^  XAA  /aV  Art  _AI*  A(9  \  ^/a9  A«f  _a9  A9  \  A9  —  /Af#  A9  _A9  a9  \A9 


*ii‘'^iii 


■‘“i-'^iv-“ii«iii>“r<“ii'^ir‘V'iv)“ii-<“o“iii-“i“ii>“ui 

^(XZ,,XZ,,-XZ,A2^^l)XZ,i,^(XZ„.Z^li-AZ,X2ii)X2,^^(XZ,X2,-X2„XZ,,)XZ„ 

■“*ll“lt-“l“lll>“u-<“o“ljr“l“il)“lll-(“t‘2r“o“u>“lv 


V«lu«a  obtAlnad  fiom  ■•van  obaarvat.iort.. 


Then,  if  k  la  A  poaltlv(  number  amallar  than  n,  ona  will  ba  able  to  avaluata 
^(n'fl)  "  aquatlona  which  ara  comprlaad  in 


'(n+l) 


-1*0  °(lc)  *  *1  *  *11 


<*0  +  *1  +  *11  * . +  *(n)>  '  *(n) 


which  muat  give  tha  aama  value  for  no  matter  which  of  tha  valuaa 

of  0,  1,  2|  3...n  ona  takaa  for  k*  If  ona  than  aolvaa  tha  aquation 


*0  *  *1  ’  *  *11  *  *111  *• 


.+  •  0 


tha  n  roota  whi^ch  will  ba  obtained  will  ba  tha  valuaa  of  p‘i  p^l  p’^i  ..tp^i 

1  ’  il’  ill  (n)' 

and  tha  valuaa  of  u^i  ate.,  of  tha  aama  aquation  will  ba  calculated 

from  tha  following  formulae*'*' i 

X  X  X  X 

(Z  -  p^^)  (Z  -  p^J^)  (Z  -  p^J)  ....  (Z  -  p^*j)  (Z  -  1) 

'^1  "  X.  X,  X.  X.  ^  X  X.  X  X 

(Pi^  -  (Pi^  -  P^J^)  (p^i  -  Pij)  ....(p^^  -  P(J))  (p^^  - 1) 


-  X  X 

**  Il[.  in  thoaa  aquatlona,  ona  dalataa  tha  factora  Z-l  and  p^  -I, 

p^^^-1,  ate.,  ona  will  have  tha  conatanta  which  ara  appropriate  to  tha 
general  term  of  tha  aariaa  dZ^,  dZ^,  dZ^^,  ate.,  i.a.,  ona  will  have  tha 
conatanta  which  ahould  multiply  ^i  .^1  ate.,  in  tha  value  of  dZ. 

Pi  .  P^i  , 
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■j 


“(Wl)  ■  ^  ■  <“l  *  "ll  +  "til  •" . ■"  “(n)> 

Obsarva  that)  aa  In  tha  pravloua  artlcla.  In  tha  coafflclanta  of  tha 

powara  of  Z  In  tha  numaratorMi  auparacrlnta  of  tha  aana  numhar  should  ba 

aubatitutadf  or  that  tha  tarms  not  containing  Z  nust  bo  multlpllad  by  Z.  and 

2  3 

thosa  containing  Z^t  ^111’  multlpllad  by  Z,  Z  ,  Z  ,  ate. 


Tha  last  valua  of  is  much  aaslar  to  calcuiato  from  tha  provlous  onoi 

which  ona  can  usa  only  for  vorlflcatlon.  If  ono  glvos  dlffaront  valuoa  to 


n,  ona  will  have  for 


Th«  abova  Is  still  In  sgrssmsnt  with  ths  fonaulss  In  my  20th  Issson 
of  analysis,  whtn  In  aach  casa  m,  and  tha  p  with  tha  hlghaat  suparacrlpt 
ara  sat  aqual  to  unity. 


Onca  tha  numbars  p^^,  *^111*  ^ 


1’  *^ll’  ‘^111’ 


ate.  ara 


obtalnad  In  this  way,  thay  can  ba  Insartad  Into  equation 


^  ■  "l  *  “it  “tl  *  “ill  “ill  “ill  * . ■“  “(n)  “(n) 


\t&2) 


which  will  satisfy  tha  2n+l  glvan  obaarvatlons,  and  can  ba  usad  to  calculata 
all  tha  valuaa  Intarmadlata  to  thoaa  observations. 


I  do  not  consider  tha  casa  where  tha  aquation 

A-  +  A .  a  +  A , .  o^  4-  . . . .+  A  •  .  a”  ■  0 

01  11  (n) 

has  4ual  or  Imaginary  roots;  In  my  20th  lesson  of  analysis,  I  have  glvan 
tha  necessary  formulas  for  finding  tha  solutions.  Xt  Is  known  that  aqual 
roots  Introduce  variable  and  rational  eoafflciants  In  the  value  of  Z,  and 
If  these  roots  ara  aqual  to  unity,  Z  will  contain  entirely  rational  terms i 
thus  tha  interpolation  formulas  which  ralata  to  rational  functions  without 
variable  dlvlslora,  or  to  parabolic  eurvas,  ara  only  a  vary  special  case  of 
those  which  I  hava  just  covered* 

I  will  now  go  on  to  applications.* 


*  Tha  applications,  while  Interesting,  ara  not  directly  relevant  to  our 
purposes  and  hava  therefore  not  bean  translated. 
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^utoQ,  J.R. ,  "SoiM  Numtrlcml  Mathoda  for  Exponantlal  Analysis  with 
Conn«ctlon  to  a  Ganaral  Idantlflcatlon  Schama  for  Llnaar  Procasaaa",  M.S. 
Thasls,  Uni varsity  of  Kantuckyt  Laxlngton,  Kantucky,  1980. 

A  ganaral  idantlflcatlon  schama  for  llnaar  processes  Is  pre¬ 
sented  and  related  to  soma  wall-known  exponential  analysis 
methods  such  as  Prony's  method  and  the  pancll-of-f unctions 
method.  An  elucidating  explanation  of  the  source  of  error 
(bias)  for  the  general  scheme  la  given.  Some  variation  of 
the  pencll-of-functlona  method  are  given  an  intense  exami¬ 
nation  with  some  revealing  numerical  examples.  A  new  ex¬ 
ponential  analysis  method*  called  the  adaptive  method*  Is 
developed  from  the  general  scheme.  Several  numerical  ex¬ 
amples  oi  the  adaptive  method  demonstrate  Its  excellent 
performance  In  the  presence  at  high  levels  of  noise.  Re¬ 
commends  t  Ions  are  made  concerning  the  use  of  the  adaptive 
method  and  new  directions  for  further  research. 


3 

4 


■j 


j 


J 


i 

i 

i 

J 
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BMCty.  L.G>,  J.D.  G«org«,  and  A.Z.  Roblnaon,  "Usa  of  tha  Coaplax  Ex** 
ponanclal  Expariaion  aa  a  Signal  Rapraaancaclon  for  Undarwatar  Acouatlc 
Callbracion»"  Journal  of  Acouatical  Soclaty  of  America.  Vol.  63,  No.  6, 
pp.  1782-1794,  Juna  19>ji. 

A  modlflad  Prony  mathod  la  davalopad  Co  provide  a  complex 
exponential  algnal  repraaencaclon  for  uae  with  underwacer 
acouatlc  calibration  waveforma.  The  approach  la  caated 
ualng  almulatad  waveforma,  and  the  almulatlon  reaulta  are 
validated  with  limited  experlmenta.  Reaulta  ahow  the  com¬ 
plex  exponential  la  an  Ideal  algnal  repreaantatlon  for 
underwater  acouatlc  calibration  bacauae  of  Ita  remarkable 
ability  to  extrapolate  calibration  waveforma  beyond  the 
actual  obaervatlon  period.  Reaulta  are  ahown  of  real  ra- 
^  clproclty  calibration  experlmenta  chat  are  valid  down  to 

23  Hz  whore  the  calibration  waveforma  are  generated  from 
a  S-aa  current  ramp  and  where  a  S-ma  obaervatlon  period 
la  uaad.  Tha  period  of  the  low-frequency  limit  la  eight 
clmea  the  obaervatlon  period  uaed  In  the  calibration. 

I  Beex,  A. A.  (Loula)  and  L.L.  Scharf,  "Rocuralve  Digital  Filter  Doalgn  Via 

Covariance  Sequence  Approximation",  IEEE  Tranaactlona  on  Acouatlea. 

Speech,  and  Signal  Proceaalng.  Vol.  ASSP-29,  No.  1,  February  IMI,  p.  51- 
57. 

Thia  paper  propoaea  a  modified  aolution  to  the  general 
^  nonlinear  leaat-aquares  approximation  problem  for  digital 

filter  deaign.  The  authors  chooae  to  approximate  tha 
covariance  aequance  associated  with  an  ideal  filter 
because  it  arises  more  naturally  in  spectnm  matching 
problems  than  dost  the  unit  ptilas  response  sequence.  The 
design  procedure  resulting  from  thair  approach  requires 
t  tha  rooting  of  an  n^^order  polynomial,  the  aolution  of 

an  nCh-order  system  of  linear  equations,  and  the  calcula¬ 
tion  of  a  (2n-l)  point  DFT  for  the  doalgn  of  an  ARMA 
(2n,  2(n-l)}  digital  filter. 

Beml,  A.J.,  "Target  Identification  by  Natural  Resonance  Estimation", 

I  IEEE  Transactions  on  Aerospace  and  Electronic  Svetema.  Vol.  AES-11,  No.  3, 

pp.  147-154,  March  1975. 

This  paper  presents  a  target  identification  method  based 
on  an  estimation  of  the  natural  fraquonclas  of  oscillation 
In  transient  radar  signatures.  The  eaphsi>is  la  placed 
^  upon  signal  modeling  and  estimation  strategy  rather  than 

relating  resonance  locations  to  physical  atructurea. 

Salient  features  of  this  Identification  method  are: 

1)  target  aspect  angle  la  not  naeded,  2)  multiple  targets 
of  the  same  type  can  be  Illuminated  simultaneously,  and 
3)  bandpass  Interrogation  pulses  can  be  uaed.  The  latter 
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faaturs  Is  compatlbla  with  existing  radar  facilities.  The 
method  Is  applied  to  soma  simulated  ti:ans£er  functional 
and  factors  affecting  estimate  accuracy  are  discussed. 

Bolllgi  6.,  and  K.J.  Langenberg,  "New  Methods  of  Nondestructive  Testing  of 
Materials  and  Their  Application  Especially  In  Nuclear  Engineering", 
Presented  at  the  International  Symposium,  laper,  Saarbrticken,  Federal 
Republic  of  Germany,  17-19  September  1979  (Text  in  German). 

A  new  method  originally  developed  for  radar  purposes  is 
discussed  for  application  to  the  acoustic  scattering 
problem.  Here  the  Helmholtz  integral  equation  Is  Laplace 
transformed  with  respect  to  time  for  the  surface  pressure 
on  a  finite  rotatlonally  synmetrlc  acoustically  hard 
scattaring  body  in  water.  The  pertinent  algebraic  equation 
is  obtained  by  Fourier  series  expansion  with  respect  to  the 
azimuthal  variable  and  by  applying  the  method  of  moments 
with  piecewise  parabolic  baae  functions  and  collocation 
testing.  The  singularities  of  the  resulting  inverse  matrix 
In  the  complex  Laplcce  plane  are  computed  for  the  sphere, 
the  oblate  and  prolate  spherlod.  These  poles  are  charac¬ 
teristic  for  the  geometry  of  the  scattering  body;,  they  can 
be  used  for  classification.  Applying  Prony's  algorithm  to 
experimental  data  similar  singularities  are  extracted. 

Brlttlngham,  J.N.,  E.K.  Miller,  and  J.L.  Willows,  "The  Derivation  of 
Simple  Poles  In  a  Transfer  Function  from  Real-Frequency  Information", 
Lawrence  Livermore  Laboratory,  UCRL-S2050,  6  April  1976. 

A  new  technique  for  extracting  simple  poles  from  real- 
frequency  transfer-function  data  Is  sunoarlzed.  The 
effects  of  varying  pole-set  parameters  are  Illustrated 
with  computed  results,  and  the  technique  Is  applied  to  the 
transfer  function  for  field  scattering  by  a  straight  wire. 

Brlttlngham,  J.N.,  E.K.  Miller,  and  J.L.  Willows,  "The  Derivation  of 
Simple  Poles  In  a  Transfer  Function  from  Real-Frequency  Information, 

Part  2:  Results  from  Real  EM  Data",  Lawrence  Livermore  Laboratory,  UCRL- 
52118,  23  August  1976. 

A  technique  for  finding  the  simple  poles  In  a  transfer 
function  F(s)  from  real- frequency  (s"ju)  data  was 
previously  reported  In  the  reference  above  and  Its 
application  demonstrated  for  a  variety  of  specific  pole 
sets.  This  procedure,  which  Is  the  analog  for  frequency 
of  Prony's  technique  for  time.  Is  further  tested  against 
specl.flc  pole  sets,  and  the  results  are  described  In 
this  report. 
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A  modification  of  tha  original  procadura,  which  uakaa 
advantage  of  tha  fact  that  ?*(»)  ■  F(8*) ,  la  praaantad. 
This  new  method  la  applied  to  electromagnetic  tranafer 
functlona.  Tha  electromagnetic  polea  thua  obtalued  are 
Invariant  with  reapect  to  excitation  and  obaervahle 
phenomena  (e.g.*  current ^  field,  etc.)  to  which  the 
tranafer  function  corraaponda,  aa  long  aa  the  syatem  is 
overdetermlned.  Varying  the  frequency  interval  over  which 
data  la  available  alao  demonatratea  the  validity  of  tha 
techxtlque. 


Buckner,  H.P.,  "Comparlaon  of  FFT  and  Prony  Algorlthma  for  Bearing 
Estimation  of  Narrow-Band  Signals  In  a  Realistic  Ocean  Environment", 
Journal  of  Acoustical  Society  of  America.  Vol.  61,  No.  3,  pp.  756-762, 
March  1977. 


Samples  of  the  quadrature  field  {Fj)  collected  on  a 
uniform  linear  array  can  be  quickly  analysed  tb  determine 
the  bearings  of  targets  by  use  of  the  FFT  or  Prony 
algorithm.  These  methods  are  compared  for  many  realisa¬ 
tions  of  a  realistic  ocean  acoustic  signal  'lua  background 
plus  nolaa  flald.  The  familiar  Monta  Ca^'lj  math' i  la 
used  to  quantify  the  algorithm  performance.  It  la  shotm 
that  bearing  accuracy  can  be  Improved  by  forming  the 
cross-sensor  field  {Gi)  where  Gi*<FjFjS4>i>.  Theoretical 
and  computational  analysis  Indlcata  that  time  averaging  of 
<FjFj*4.1>  not  only  reduces  the  effects  of  noise,  but  also 
reduces  Interactions  of  i/avas  In  the  signal  and  background 
fields  resulting  In  better  array  performanca. 


I 


I 


I 
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Carr,  John  W. ,  "An  Analytic  Investigation  of  Transient  Synthesis  by 
Exponentials",  Thesis,  Massachusetts  Institute  of  Technology,  1949. 

This  paper  Investigates  least-mean- square  approximations 
to  transient  functions  In  the  time  domain,  offering  a 
solution  to  the  problem  of  synthesis  of  networks  to  give 
specified  transient  responses.  Various  exact  solutions 
to  this  problem  are  Investigated,  none  of  which  yields  a 
solution  using  a  high-speed  electronic  digital  computer. 
Further  Investigation  is  mde  of  useful  "approximations 
to  the  best  approximation."  Analysis  begins  with  the 
Prony  method  of  undetermined  coefficients  for  ordinary 
linear  differential  equations  with  constant  coefficients, 
extending  this  method  to  the  case  of  Integral  equations 
as  well.  This  enables  functions  without  continuous 
derivatives  to  be  approximated  by  a  set  of  decreasing 
exponentials,  either  real  or  complex. 

Chan,  L.C.,  and  L.  Peters,  Jr.,  "Subsurface  Electromagnetic  Mine 
Detection  and  Identification",  Ohio  State  University  ElectroSclence 
T.aboratory,  Technical  Report  784722-1,  November  1978. 

The  goal  of  this  research  Is  to  adapt  a  conmercially 
available  radar  system  developed  for  the  detection  of 
pipes  to  the  detection  of  shallow  objects  and  further  to 
establish  a  target  Identification  capability.  The  target 
Idontlficatlon  technique  selected  for  this  purpose  (from 
among  several  being  pursued  at  the  ElectroSclence 
Laboratory)  Is  designated  the  predictor-correlator 
scheme.  This  Identification  scheme  has  already  been 
extremely  successful  In  that  a  "mlnellks"  target  has  been 
separated  from  a  set  of  false  targets.  The  best  designed 
and  completely  tested  system  to  date  yields  lOOZ  Identifi¬ 
cation  of  the  "mlnellke"  target  and  a  zero  false  alarm 
rate  of  any  of  the  false  targets.  There  remains  substan¬ 
tial  Improvements  (some  of  which  are  nearlug  completion) 
that  can  be  Incorporated  in  both  the  radar  system  and  In 
the  target  Identification  techniques.  However,  the  target 
Identification  capability  discussed  In  this  report 
represents  such  a  substantial  advance  that  It  is  being 
published  in  Its  present  form. 

Chan,  L.C.,  "Subsurface  Electromagnetic  Target  Characterization  and 
Identification",  Ohio  State  University  ElectroSclence  Laboratory, 
Technical  Report  784722-3,  April  1979. 

A  method  for  subsurface  radar  target  characterization  and 
identification  Is  described.  This  method  characterizes 
subsurface  radar  targets  by  their  complex  natural 
resonances  which  are  extracted  directly  from  their 


backscattered  time-domain  waveforms.  The  difference 
equation  coefficients  associated  with  the  complex 
resonances  are  then  used  In  the  predictor-correlator  for 
target  Identification.  Both  the  characterization  and 
Identification  processes  are  extensively  tested  with  real 
radar  measurements  and  found  ro  yield  practical  target 
Identification  performance.  The  target  Identification 
process  Is  simple  and  Involves  only  simple  algebraic 
operations.  Based  on  the  Identification  process,  a  "first- 
generation"  microcomputer  Identification  radar  system  Is 
Implemented  for  target  Identification  In  real  time.  This 
radar  system  la  found  to  yield  practical  Identification 
performance. 

Chan,  L.C.,  D.L.  Moffatt,  and  L.  Peters,  Jr.,  "A  Characterization  of 
Subsurface  Radar  Targets",  Special  Issue  on  Exploration  Geophysics, 
Prop,  of  IEEE.  July  1979. 

The  capability  of  subsurface  target  identification  at 
shallow  depths  has  been  demonstrated  using  an  electro¬ 
magnetic  video  or  baseband  pulse  radar.  Real  radar 
measurements  were  collected  for  five  targets  at  a  depth 
of  5  cm  (2  in)  in  various  ground  conditions.  These 
measurements  were  processed  for  target  characterization 
and  Identification.  Identification  performance  based  on 
a  single  radar  observation  was  evaluated.  The  Identifi¬ 
cation  process  requires  only  simple  algebraic  operations 
and  thus  offers  the  potential  of  real-time  on-locatlon 
Identification  of  subsurface  targets.  Pole  patterns  are 
extracted  from  the  backscattered  waveforms  using  Prony's 
method. 

Cho,  K.S.  and  J.T.  Cordaro,  "Calculations  of  the  Singularity  Expansion 
Method  Parameters  from  the  Transient  Response  of  a  Thin  Wire",  AFWL 
Interaction  Note  379,  September  1979. 

The  problem  of  determining  the  singularity  expansion 
method  (SEM)  parameters  from  transient  thln-wlre  data 
is  examined.  A  computer  code  Is  used  to  generate  re¬ 
sponse  data.  The  SEM  parameters  are  computed  from  these 
data  using  Prony's  method.  For  noisy  dat;^,  it  is  shown 
that  the  parameter  values  can  be  improved  by  signal 
averaging. 

Chuang,  C.W.  and  D.L.  Moffatt,  "Complex  Natural  Resonances  of  Radar 
Targets  Via  Prony's  Method",  Ohio  State  University  ElectroSclence 
Laboratory,  Technical  Report  3424-3,  April  1975. 
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It  haa  baen  shown  that  axtractlon  of  complex  natural 
resonances  and  residues  from  exact  transient  responses  via 
Prony's  method  Is  both  efficient  and  accurate.  However, 
in  the  art  of  radar  Interrogation  the  exact  transient 
response  is  difficult  to  obtain  and  we  must  be  content  with 
an  approximate  one.  Knowing  the  complex  natural  resonances 
of  objects  Is  essential  to  one  type  of  radar  target 
discrimination.  In  this  paper  we  consider  the  use  of 
Prony's  method  In  extracting  complex  natural  resonances 
from  approximate  backscattered  ramp  responses  of  radar 
targets . 

A  ramp  waveform  la  used  In  radar  target  discrimination 
because  a  good  approximation  to  the  time  domain  response 
can  be  obtained  with  only  a  few  harmonically  related 
frequency  domain  responses  via  Fourier  synthesis.  By 
applying  Prony's  method  to  this  approximate  ramp  response, 
a  few  dominant  resonances  which  are  close  to  the  origin  of 
the  complex  frequency  plane  can  be  extracted.  This  Is 
quite  encouraging  because  for  the  purpose  of  radar  target 
discrimination  only  a  few  dominant  modes  are  needed. 

Objects  considered  are  prolate  spheroids,  simple  wire 
aircraft  models  and  some  realistic  aircraft  models.  The 
accuracy  of  this  method  Is  found  to  be  good  when  compared 
to  the  results  obtained  by  other  techniques. 

Chuang,  C.W.  and  D.L.  Moffatt,  "Hatural  Resonances  of  Radar  Targets  Via 

Prony's  Method  and  Target  Discrimination",  IEEE  Transactions  on 

Aerospace  and  Electronics  Systems,  Vol.  AES**  12,  Mo.  5,  pp.  583**S89, 

September  1976. 

The  dominant  complex  natural  resonances  of  radar  targets 
are  obtained  via  method  applied  to  calculated  and 

measured  backscatter^c  ra)?(p  response  waveforms.  Subject 
targets  are  spheres,  simple  wire  models  of  straight  and 
swept  wing  aircraft,  end  realistic  models  of  modern  fighter 
aircraft.  It  Is  demonstvaced  ilit  when  the  backscattered 
ramp  response  waveforms  are  obtained  via  Fourier  synthesis 
of  limited  spectral  range  harmonic  scattering  data,  some 
resonance  locations  at  variance  with  those  obtained  from 
reaction  integral  equation  search  procedures  are  obtained. 

It  Is  also  shown,  however,  that  the  Prony  deduced  resonances 
can  be  used  successfully  in  predlctor**correlator  target  dis¬ 
crimination. 


Cordaro,  J.T.,  "A  Note  on  Representing  e  Transient  Waveform  by  a  Finite 
Sum  of  Complex  Exponentials",  AFWL-TR‘*77-46,  Air  Force  Weapons  Laboratory, 
July  1977.  Also  Air  Force  Weapons  Laboratory  Mathematics  Notes  #46. 

The  problem  of  representing  an  observed  transient  waveform 
by  a  finite  sum  of  complex  exponentials  is  considered.  A 
least  squares  technique  that  involves  overfitting  the  model 
signal  and  then  eatlmatlng  the  correct  order  is  applied  to 
simulated  data  as  well  as  data  recorded  on  photographs.  In 
all  the  examples,  the  fitted  waveform  approximates  the 
original  waveform  well  when  both  functions  are  plotted  on 
the  same  graph. 

Cordaro,  J.T.,  "Comparison  of  Three  Techniques  for  Calculating  Poles  and 
Reaiduaa  from  Experimental  Data,"  AFWL  Mathematics  Notes  #61,  August  1978. 

Three  pole  and  residue  calculation  techniques  are  compared 
by  applying  the  methods  to  simulated  data.  Each  method 
works  wall  when  the  data  consist  of  a  damped  sinusoid  plus 
noise.  But  when  the  data  include  effects  from  pulser 
asynchronlam  and  ground  reflection,  accurate  pole  calcu¬ 
lation  is  more  difficult.  For  the  methods  to  produce 
accurate  results  under  some  EMP  test  conditions,  it  is 
necessary  to  use  incident  field  data  and  information  about 
the  SEM  coupling  coefficients  and  natural  modes  of  the 
body  exposed  to  EMP. 

Cordaro,  J.T.,  "Pole  Measurements  for  the  ATHAMAS  Pipe  Test",  AFWL 
Mathematics  Motes  #56,  August  1977. 

In  the  ATHAMAS  pipe  test,  surface  current  and  charge  density 
response  data  were  recorded  for  a  metal  cylinder  illuminated 
by  EMP.  Poles  and  residues  describing  some  of  these  data 
have  been  computed  using  an  iterative  technique.  Some  poles 
can  be  Identified  as  corresponding  to  natural  moden  of  the 
cylinder.  Others  match  the  Incident  field  poles.  Most 
poles  in  the  incident  field  are  shown  to  be  due  to  ground 
reflection. 


d«  Flgutlrtdo,  Rul  J.P.  and  C.L.  Hu,  "Appllcaclon  of  a  Fraquancy  Domain 
Prony  Hathod  co  Wlda  Bandwidth  Radar  Slgnatura  Claaaificacion,"  Daparcmant 
of  Elactrlcal  Enginaarlng,  Rica  Unlvarilty,  Houaton,  Taxait  Tachnlcal 
Report  EE  7908,  Saptambar  1979. 

Lat  X  and  y  ba  slgnala  (l.a.,  raal-valuad  functlona  of  tlma) 
of  finite  duration  and  anarjy.  Thli  paper  davalopa  a  fra- 
quancy  domain  Prony  approach  for  Interpolating,  or,  in 

general,  approximating  y(t)  by  a^  x(t-T^)S  where 

a  ■  (a].,  ...  ,  aj()  and  X"  (tj^,  ...  ,  tj^)  are  real 

M-vector  parametera.  If  x  and  y  denote  the  Fourier 
traniforms  of  x  and  y,  then  H((d)  ■  Y(u)/X(u)  la  a  linear 
combination  of  exponential ev^h  coafflclanti  ai;  and 
exponente  -jTiw,  whore  j  ■  / -1  .  The  Prony  method  la 
uaod  to  determine  ai  and  ti  from  the  aamplea  of  H.  Theao 
aamplea  are  obtained  trom  the  (time-domain)  aamplea  of  x 
and  y  by  meana  of  the  OFT.  Both  the  caaea  In  which  theao 
aamplea  are  nolae-freo  and  contaminated  by  nolae  are  con- 
aldered.  Special  modlf Icatlona  of  the  conventional  Prony 
procedure  to  take  Into  account  the  fitting  of  complex  data 
are  preaentad.  Theao  reaulta  are  then  extended  to  the 
caae  where  y  reaulta  from  a  linear  combination  of  aeveral 
known  algnala  xi,  each  differently  delayed.  The  caae  In 
which  y  reaulta  from  the  convolution  of  a  known  aignal  x 
with  an  unknown  Impulao  raaponae  h  la  alao  conaldered. 

A  connection  between  Wiener 'a  Tauberlan  thaorema  and  the 
frequency  domain  Prony  approximation  are  preaentad  here. 

de  Flgualredo,  Rul  J.P.  and  C.L.  Hu,  "A  Fourler-Prony  Tauberlan  Approach 
tu  che  Analyala  of  a  Mixture  of  Delayed  Slgnala,"  Department  of  Elactrlcal 
Engineering,  Rice  Unlvaraity,  Houaton,  Texaa,  Technical  Report  EE  7907, 
September  1979. 

A  frequency  domain  Prony  approach  la  preaentad  for  ex¬ 
tracting  featurea  of  return  algnala  from  targeta  illumi¬ 
nated  by  wide  bandwidth  (abort  pulae)  radar.  Theoretical 
detalla  pertaining  to  thla  approach  are  daacrlbed  in  the 
above  rnfarenca.  The  featurea  conalat  of  the  relative 
delaya  and  reflection  coefflclenta  pertaining  to  acatterlng 
centera  on  the  target  repreaentlng  differently  ahaped  re- 
glona  on  the  target  aurface.  The  dlmenelonallty  of  the 
feature  vectora  thua  ccnatructad  la  very  low  (leaa  than  ten). 
Moreover,  when  uaed  in  the  claaalflcation  of  targeta  by  a 
neareat  neighbor  claaalflcation  atratagy,  auch  feature 
vectora  permit  accurate  dlacrlmlnatlon  between  targeta  that 
do  not  differ  much  in  ehape;  and  alao  they  are  in  a  large 
meaaure  Inaenaltlve  to  nolae.  The  reaulta  preaentad  were 
corroborated  by  computer  aimulatlona  performed  on  the  data 
baae  created  by  the  coherent  X-band  ahort  pulae  (0.5  nano¬ 
second)  radar  at  the  Fort  Worth  operation  Radar  Range  of 
General  Uynomlca  Convalr  Aeroapace  Diviaion. 
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Dudl«y,  D.G.,  "Fitt:ing  Noisy  Dscs  with  s  Complsx  Expontnclsl  S«rls»," 
UwTsnct  Livsrmors  Lsborseory»  UCRL-S2242,  7  Msrch  1977.  Also  AFWL 
Mschsmseics  Notts  //SI. 

Tht  procsss  of  fitting  noisy  dtts  with  s  complsx  tx- 
pontntltl  strlts  Is  consldtrtd.  An  sntlysls  of  tht 
cltsslctl  Prony  mtthod  shows  thtt  ttttmpts  to  obtain 
accuratt  anawars  In  tht  prtstnct  of  nolaa  matt  with 
consldtrablt  dlfflcultltt.  Rtcastlng  tht  probltm  as 
ont  of  systam  Idtntlflcatlon  and  applying  tht  Itast- 
squarta  mathod  Itada  to  blaatd  tstlmacts  of  tht  para- 
matara  In  tha  charaetarlstlc  aquation  for  tha  systam. 

Tha  difficulty  cantara  around  tha  fact  that  nalthar 
classical  Prony  nor  least  aquaras  Involvas  any  analysis, 
Idtntlflcatlon,  smoothing,  or  filtering  of  tha  nolaa. 

Soma  tachnlquas  for  Improving  results  art  discussed, 
and  an  Itaratlva  ganarallsad  laast-squaraa  procadura, 
which  leads  to  a  nolaa  filter,  la  racommandad  for  further 
study. 

Dudley,  D.G. ,  "Simple  Discrete  System  Modelling  In  Hilbert  Space,"  Uni¬ 
versity  of  Arisons,  1  February  1978. 

Tha  design  of  a  simple  discrete  systam  model  for  a 
aaquanca  of  data  is  given.  Tha  modal  la  shown  to  be 
compatible  with  approximation  In  Hilbert  space.  Tha 
minimum  norm  solution  Is  obtained.  Tha  result  la  a 
structure  which  accommodates  many  of  tha  results  of 
system  Identification  and  linear  prediction.  The 
effects  of  error  are  briefly  mentioned. 

Dudley,  D.G.,  "Parametric  Modeling  of  Transient  Electromagnetic  Systems," 
^2jl^o^„jclence,  Vol.  14,  No.  3,  pp.  387-396,  May-June  1979. 

The  modeling  of  a  transient  electromagnetic  system  as 
a  single-input,  single-output,  linear,  casual  process 
Is  considered.  Two  canonical  forms  era  discussed t  one 
for  impulse  function  Input  and  tha  other  for  more  general 
Input.  Expressions  are  derived  for  errors  In  the  em  l- 
matlon  of  the  model  parsmaters  caused  by  errorc  In  the 
process  data.  It  is  shown  that  Prony 's  method  Is  a 
special  case  of  the  Impulse  function  modal  and  that 
serious  errors  occur  in  attempts  to  apply  tha  method  in 
a  noisy  environment.  It  Is  also  shown  that  least  squares 
estimation  of  the  modal  peramatars  produces  blasad  re¬ 
sults  and  that  attempts  to  overcome  this  problem  have 
led  to  nonlinear  equation  sets.  Some  recommendations 
are  made  concerning  noise  Identification  and  filtering 
in  order  to  improve  results. 
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Duh*m«l.  P.  tnd  J.C.  R«ult,  '‘AutoMitlc  T.it  C«Mrttlon  TtchniquM  for 
An»iog  Circuit*  «nd  Systmii  A  R*viiw,”  IBBB  9n  <;U9Ut£l 

■nd  Sviit«m*.  Vol.  CA8-26,  No.  7,  pp.  411-440,  July  1979. 

Thli  p*p«r  t*  both  a  raviaw  and  an  aaaaaamant  ot  tach- 
nlquaa  praaancly  avallabla  for  automatic  taat  ganaration 
for  analog  ayatama  tachniquaa  for  catagoriaatlng  tha 
problama  of  automatic  taating  (daflnitiona,  fault*  In 
analog  ayatama,  diffaranc  typaa  of  taata, 
and  dlagnoala  procaduraa),  charactarlaatlon 
modaa  of  analog  ayatama,  and  tha  main 

of  automatic  taat  aquipmant  can  b*  propoaad.  Savaral  tach- 
niquaa,  raapactivaly,  procaading  from  approachaa  baaad  on 
datarminiatio  and  probahiliatic  aatimation,  taxonomlcal 
and  topological  analyaaii  ara  datailad.  Tachniquaa  apaciflo 
to  llnaar  ayatama  (aavaral  of  thorn  balonging  to  tha  abova 
thraa  catagoriaa)  ara  daalt  with  in  a  aaparata 
Prony'a  mathod  la  daacribad  aa  a  mathod  of  tranafar  function 
datarmination. 
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Fl«tch«r,  D.L.  *nd  C.N.  Wtygtndt,  "A  Dlgitui  Method  of  Trtnifor  Function 
r.*lcuX*tion,"  lEBE  Tr>ni.  Clr.  Thtory.  Vol.  CT-18,  pp.  185-187,  January  1971. 

A  maehod  li  davalopad  that  allowa  tha  dlract  calculation 
of  tha  laaat-iquaraa  aatimataa  of  tha  polaa  and  raaiduaa 
of  tha  tranafar  function  of  a  llnaar  lyatam  whan  tha  input 
and  output  of  tha  ayatam  la  rapraaantad  In  a  aaquantlal 
nutaarical  format.  Tha  mathod  la  particularly  adaptable  to 
Implanantation  on  a  digital  computer  and  la  qulta  efficient 
for  problama  of  modarata  alaa.  Tha  method  praaantad  can  be 
conaldarad  an  axtanalon  to  Prony'a  method  aa  applied  to 
tranafar  function  uynthaala. 

Frooma,  D.A.,  V.A.J.  van  Lint,  T..  Stattnar,  R.W.  Maegurn,  and  O.M. 

Goldataln,  "Explodlng-Wlra  Photon  Taatlng  of  Skynat  Satellite, "  IEEE 

Vol.  NS-23,  No.  6,  pp.  1349-135T; 


Thla  paper  praaanta  a  comparlaon  of  part  of  tha  data  ra- 
aultlng  from  a  aariaa  of  photon,  alactrlcal  and  caloulatlonal 
axarclaaa  rm  a  atructural  modal  of  tha  Skynat  1.  Calculationa 
made  for  comparlaon  with  tha  teat  data  ware  of  two  format  coda 
aimulationa  and  analytic  modala.  Tha  code  calculationa  wara 
tha  moat  detailed  and  Included  rapraaantatlona  of  tha  major 
thraa-dlmanaional  gaomatrlcal  faaturaa  that  wara  apparent  on 
tha  Skynat  Structure  Modal  <S8M).  For  the  muat  part,  pra- 
dlctlona  wara  made  at  high  flux,  uaing  a  praacrlbad  aourca 
and  a  generic  pulaa.  Analytic  modala  wara  uaad  to  under atand 
tha  raaulta  produced  by  tha  coda  calculationa  and  tha  axparl- 
manta.  Simple  analytic  circuit  modala  rapraaantlng  alactro- 
magnatlc  modaa  of  tha  ayatam  ahow  promiaa  of  being  raaaonably 
accurate  pradlotiva  maana  for  rapid  computation  of  tha 
alaotromagnatlc  raapnnaa. 


Gardner,  C.D. ,  "Expali-Fortran  Coda  for  Expontnclal  Approximation  by 
LitiC  Squacai,"  Lawttnea  Radiation  Laboratory,  UCRL-14S41,  Rev.  1, 
Unlvarilty  of  California,  Llvtrmora,  March  1970. 

EXPALS  Is  a  FORTRAN  program  for  obtaining  a  wslghtsd  Isast- 
squarts  fit  of  a  sum  of  axponantials  to  bt  run  on  CDC  6600 
or  IBM  7094  computers.  This  routine  fits  a  linear  com¬ 
bination  of  real  exponential  decay  functions  by  the  Prony- 
Householder  iteration  method.  This  report  presents  the 
mathematics  and  progroioming  necessary  to  run  EXPALS.  Also 
included  ere  output  of  examples  to  test  the  program. 

Gavel,  D.T.,  J.V.  Candy,  and  D.L.  Lager,  ’’Parameter  Estimation  from  Noisy 
Transient  Electromagnetic  Measurement,”  Lawrence  Livermore  Laboratory, 
Submitted  for  Nuclear  EKP  Meeting,  Anaheim,  California,  5-7  August  1980. 

This  paper  describes  the  performance  of  three  different  signal 
processing  techniques  applied  to  parsmeteriss  a  body  from 
noise  experimental  electromagnetic  transient  response  date. 

The  techniques  range  from  the  well-known  Prony  method  to  the 
more  sophisticated  extended  Kalman  filter  end  finally  to  the 
highly  sophisticated  maximum  likelihood  identifier.  Comparison 
is  made  of  the  performance  of  these  algorithms  end  a  discussion 
of  their  tradeoffs  is  given. 

Gerck,  E. ,  "Prony  Method  Revised,”  j^22li£^jO££jjC£,  Vol.  18,  No.  18, 
pp.  3075,  15  September  1979. 

The  Prony  method  described  in  its  application  to  optical 
resonators  by  Siegmeit  and  Miller  has  since  then  deserved 
some  attention.  This  letter  shows  that  the  Prony  method 
for  symmetric  matrix  is  ell  that  is  necessary  for  any 
scaler  resonator  calculation,  regardless  of  misalignments 
or  misfiguring  of  mirrors.  This  is  accomplished  by  simple 
modifications  of  the  original  algorithm.  In  this  way  s 
simpler  end  faster  procedure  is  presented  only  by  con¬ 
sidering  the  correct  scaler  product  to  be  used.  The  pro¬ 
posed  method,  however,  could  not  possibly  apply  directly 
in  the  cases  for  which  the  polarisation  of  the  beam  is 
changed  aa  in  "optical  diodes”  or  sxicons. 

Goodman,  D.K. ,  "Determination  of  the  Number  of  Terms  Necessary  for  a 
Class  of  Approximation  Procedures,"  in  Data  Representation.  Edited  by 
H.S.  Andersen  end  M.R.  Osborne,  University  of  Queensland  Press,  pp. 77-93, 
1970. 


This  presentation  is  concerned  with  the  following  areas t 
(1)  the  order  determination  task  when  the  basis  functions 
are  exponentials;  (2)  a  tutorial  Introduction  to  a  new 
arithmetic  designed  for  digital  computation  which  pro¬ 
duces  automatic  error  bounding;  end  (3)  a  solution,  both 
theoretical  end  practical,  to  the  problem  in  (1)  end  the 
exponential  fitting  procedure. 
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H«nd«rson,  T.L,,  "Matrix  Mathoda  for  Datamlnlng  Syatam  Polaa  from 
Tranalant  Raaponaa",  Unlvaraity  of  Kantucky,  Elactromagnatlca  Raaaarch 
Raport  80-1,  May  1980,  alio  AFVfL  Mathamatlca  Nota  //66. 

Tha  problam  traatad  la  that  of  Idantifying  tha  polaa  of  a 
flnita  ordar  ayatam  by  obaarvlng  ita  tranalant  dacay  aftar 
caaaation  of  input,  for  a  llmltad  tima,  uaing  (poaaibly) 
fflultipla  obaarvatlon  polnta  and  axparimantal  rapatitlon. 

Varioua  approachaa  ara  atudlad,  having  tha  common  charac- 
tariatic  that  a  homoganaoua  matrix  aquation  muat  ba  aolvad. 

Savaral  tachniquaa  that  hava  baan  givan  acant  attantion  in 
tha  litaratura  ara  conaolidatad  into  tha  traatmant,  togathar 
with  naw  raaulta  including  an  analytical  traatmant  of  tha 
conaaquancaa  of  aaauming  an  axcaaalvaly  high  ayatam  ordar, 
darivation  of  a  atatiatically  unbiaaad  aatimata  for  an 
intarmadiata  paramatar  in  tha  aolution,  naw  thaorama  on 
arror  af facta,  a  racipa  for  affactiva  uaa  of  tha  aingular 
valua  dacompoaition,  a  naw  mathod  for  auppraaaion  of 
axtranaoua  polaa,  an  alucldating  darivation  and  axtanaion  of 
tha  mathod  of  Jain,  a  naw  form  of  tha  problam  wharain 
tha  ayatam  polaa  ara  aiganvaluaa,  and  a  atudy  of  tha 
ralationahip  batwaan  varioua  pola  idantif Ication  mathoda. 

Handaracn,  T.L.,  "Oaomatric  Mathoda  for  Datarmlning  Syatam  Polaa  from 
Tranalant  Raaponaa",  to  ba  publiahad  by  ISES  Tranaactlon  on  ASSP. 

Tha  problam  traatad  ia  that  of  idantifying  tha  polaa  of  a 
finlta  ordar  ayatam  by  obaarvlng  ita  tranalant  dacay  aftar 
caaaation  of  input,  for  a  llmltad  tlma,  uaing  (poaaibly) 
multlpla  obaarvatlon  polnta  and  axparimantal  rapatitlon. 

Varioua  approachaa  ara  atudlad,  having  tha  common 
charactarlatlc  that  a  homoganoua  matrix  aquation  muat  ba 
aolvad.  Savaral  tachniquaa  that  hava  baan  givan  acant 
attantion  In  tha  litaratura  ara  conaolidatad  by  a 
gaomatric  traatmant,  togathar  with  naw  raaulta  including 
an  analytical  traatmant  of  tha  conaaquancaa  of  aaauming  an 
axcaaalvaly  high  ayatam  ordar,  darivation  of  a  atatiatically 
unbiaaad  aatimata  for  an  intarmadiata  paramatar  in  tha 
aolution,  naw  thaorama  on  arror  af facta,  a  racaipa  for 
affactiva  uaa  of  tha  aiagular  value  dacompoaition,  and  a 
naw  mathod  for  auppr'-aolon  of  axtranaoua  polaa. 

Hillar,  J.,  "Exponantial  Fitting  and  Ganaral  Tachniquaa  of  tha  Prony 
Form",  in  Data  Rapraaantation.  adltad  by  Andaraan  and  Oaborna,  Unlvaraity 
of  Quaanaland  Praaa,  pp.  62-76,  1970. 

Tha  intantion  of  tha  firat  part  of  th.  papar  ia  to  poaa  a 
aarlaa  of  quariaa  and  to  praaant  aoma  axamplaa  that  aarva 
to  ahow  tha  ralavanca  of  axponantlal  fitting.  It  ia 
damonatratad  that  tha  juatlfication  for  axponantlal 


fitting  Is  oftsn  poor.  Nonsthslsas  it  is  not  always  assy 
to  posa  tha  "corract"  fitting  problam  and  axponantial 
fitting  may  ba  raasonably  invokad  by  dafault.  Tha  sacond 
part  of  tha  papar  diacuaaas  tha  Prony  procadura  in  tanas 
of  commuting  oparators.  Tha  aim  of  tha  traatmant  hare  is 
to  point  out  tha  quastiona  that  should  ba  askad  so  that 
ganaral  oparator  claasaa  may  ba  handlad. 

Holt,  John  N.  and  Robyn  J.  Antill,  "Datarmining  tha  Numbar  of  Tanas  in  a 
Prony  Algorithm  Exponantial  Fit,"  Mathamatieal  Bloselancas.  36,  pp.  319> 
332,  1977. 

Tha  problam  of  datarmining  tha  corract  numbar  of  tanas  to 
taka  In  an  axponantial  sum  approximation  to  squally  spacad 
data  is  considarad.  Tha  Prony  algorithm  is  usad,  togathar 
with  singular  valua  dacompoaition  tachniquas.  A  critarion 
la  aatabliahad  which  is  uaaful  in  aatimating  tha  numbar  of 
tarma  whan  tha  data  ara  sub j act  to  normally  distributed, 
zaro  maan  noiaa. 

Houaaholdar,  A.S.,  "On  Prony'a  Mathod  of  Fitting  Exponantial  Dacay  Curvaa 
and  Multipla-Hlt  Survival  Curvaa,"  Oak  Ridga  National  Laboratory,  ORNL- 
4SS,  Oak  Ridga,  Tannassaa,  Fabruary  19S0. 

This  papar  first  daacribas  Prony 's  mathod.  Tharaaftar  it 
is  shown  how  ona  can  procaed  to  obtain  a  valid  laaat 
squaras  fit,  and  finally  a  critarion  for  choosing  tha 
number  of  axponantials  raquirad  is  given.  This  is  tha 
source  of  tha  Householder  orthogonalization  procadura. 

Hudson,  H.G.  and  D.L.  Lager,  "Observations  on  tha  Operation  of  tha  SEMPEX 
Coda,"  Lawranca  Livermore  Laboratory,  UCID-17440,  22  Saptambar  1976. 

The  authors'  observations  on  cha  use  of  tha  SEMPEX  computer 
coda  ara  given.  The  issuas  of  proper  sampling  rate  and 
total  aampling  interval  are  addressed.  The  effect  that 
noisa  in  tha  data  .'s  on  tha  results  and  cha  ability  of 
autocorrelation,  truncation  filtering,  and  cha  sliding- 
window  technique  to  reduce  chase  affacta  are  also  discussed. 


Jain,  V.K. ,  "Raprasantatlon  of  Saouancaa.'*  IEEE  Transactions  on  Audio  and 
Elsctroacoustlcs.  Vol.  AU-19,  No.  3,  pp.  208-215,  Saptsmbsr  1^71.  ~ 

In  this  papsr,  tha  spaca  of  squara  suimnabla  raal  saquancas 
Is  consldarad.  Insight  Is  davalopod  Into  tha  structura  of 
tha  spaca  and  a  mathod  for  raprasantatlon  of  Its  alamants 
Is  dascrlbad.  Tha  mathod  la  shown  sultabla  for  on-llna 
Implamantatlon  and  la  nonitaratlva;  It  can  ba  usad  for 
wavaform  daslgn,  signal  faacura  axtractlon,  and  dlscrata- 
tlma  systam  Idantlflcatlon.  Tha  applications  ara  axampllflad 
by  thraa  simulation  studias. 

Jain,  Vljay  K.,  Tapan  K.  Sarkar,  and  Donald  D.  Walnar,  "Extanslon  of 
Pancll-of-functlona  Mathod  to  Raversa-Tima  Procaaslng  with' First-Order 
Digital  Flltars,"  Dapartmsnt  of  Elactrlcal  Englnaarlng,  Rochastar 
Instltuta  of  Tachnology,  Rochestar,  Naw  York,  Tachnic^l  Raport  No.  3, 
August  1980. 

In  this  prasantatlon,  tha  data  signal  Is  procassed  In 
ravarsa-tima  by  a  cascada  of  first  ordar  digital  flltars 
to  ylald  a  family  of  Information  signals.  Tha  Gram  matrix 
of  thasa  information  signals  Is  shown  to  contain  tha 
assantlal  Informatlo.i  on  tha  polas  of  tha  signal.  Tha 
antlra  procadurc  of  tha  application  of  pancll-of-functlon 
mathod  Is  thus  nonitaratlva.  Examplas  prasantad  damon- 
strata  (1)  nolsaworthinasa  In  tha  raprasantatlon  problam 
whan  data  ara  corrupted  by  noise,  and  (11)  tha  affactivenass 
of  tha  mathod  In  the  approximation  problam. 
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KanB&l«r,  D.W.,  "Elxlsccnce  of  B«Bt  Approximations  by  Sums  of  Exponantials," 
Journal  of  Approximation  Thsorv.  Vol.  9,  No.  1,  pp.  78-90,  Septembar  1973. 

This  thssis  studiss  tha  problem  of  approximation  by  sums 
of  exponentials  and  contains  a  discussion  of  Prony's 
method . 

Kammler,  D.W.,  "Prony's  Method  for  Completely  Monotonic  Functions," 

Journal  of  Mathematical  Analysis  and  Applications.  Vol.  57,  No.  3, 
pp.  560-570,  1  March  1977. 

Prony's  method  for  constructing  an  approximating  exponential 
sum  requires  one  to  select  2n  uniformly  spaced  points,  solve 
a  difference  equation  to  find  ci,  ...  ,  Cn,  find  the  roots 
of  a  polynominal  and  take  a  log  to  find  the  exponents 
111  •••  I  In*  procedure  cannot  be  used  to  produce  an 

exponential  sum  for  an  arbitrary  FeC[0,»l  but  then  one 
cannot  reasonably  expect  to  model  an  arbitrary  continuous 
function  with  such  an  exponential  sum.  This  paper  proceeds 
under  the  additional  hypothesis  that  F  is  completely  nono¬ 
tonic  on  [0,<*].  It  is  shown  that  Prony's  method  always 
yields  an  exponential  sum  which  for  sufficiently  large  n 
can  be  made  arbitrarily  cose  to  F  (in  the  uniform  norm)  if 
and  only  if  F  is  completely  monotonlc  on 

Kelso,  J.M.,  "Measuring  the  Vertical  Angles  of  Arrival  of  HF  Skywave 
Signals  with  Multiple  Modes,"  Radio  Science.  Vol.  7,  No,  2,  pp,  245-250, 
February  1972. 

The  measurement  of  the  vertical  angle  of  arrival  of  HF 
skywave,  multimode  signals  generally  requires  either 
highly  complicated  mathematical  manipulations  or  observing 
tachnlques  Chat  necessitate  a  waiting  period  until  the 
relative  phases  among  the  modes  either  reach  a  desired  con¬ 
dition  or  pass  through  a  desired  range  of  variations.  This 
paper  utilizes  an  extension  of  a  very  old  method  (the  Prony 
method,  c.  1795)  for  approximating  sums  of  exponential  func¬ 
tions.  The  technique  uses  relatively  simple  calculations  that 
are  based  on  ^an  instantaneous  (or  nearly  so)  set  of  samples 
from  a  uniformly  spaced,  vertically  disposed  array  of  antenna 
elements.  Unlike  many  other  methods,  this  technique  requires 
neither  time  variations  of  phase  nor  a  perfectly  reflecting 
ground  plane,  and  is  valid  for  arbitrary  polarization  of  the 
(similar)  antenna  elements  and  for  any  general  wave  polariza¬ 
tion.  Experimental  results  from  some  field  tests  of  the  pro¬ 
cess  for  vertical  angle  esclmatlon  are  presented. 

Khatwanl,  K.J.,  and  J.S.  Bajwa,  "Identification  of  Linear  Time- Invar lane 
Systems  Using  Exponential  Signals,"  IEEE  Transactions  on  Automatic  Con¬ 
trol.  pp.  146-148,  February  1975.  ——————————— 
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A  new  method  o£  system  identification  for  linear  time 
invariant  systems  based  on  Prony's  method  of  exponential 
interpolation  is  suggested  in  this  correspondence.  This 
method  of  identification  uses  an  exponential  signal  for 
system  excitation  and  has  the  advantages  of  simplicity, 
independent  determinatlou  of  numerator  and  denominator 
of  the  system  transfer  without  any  difficulty.  The 
simplicity  of  the  method  has  resulted  because  of  the 
simple  but  particularly  suitable  formulation  in  the  state 
variable  form  used. 

Kramer,  H.  and  A.  Lacroix,  "Praktische  Beref;hnung  von  Exponentlal- 
approxlmatlonen  mit  dam  Prony-Verf ahren" ,  Anaewandte  Infortaatik.  Vol.  3, 
pp.  121-126,  1972. 

Prony's  method  is  described  for  the  approximation  of  discrete 
sets  of  points  by  exponentisl  functions  with  complex  argument.' 

The  determination  of  the  degree  of  the  approximating  function 
is  discussed  in  detail.  It  is  shown  by  examples  how  Input 
•^.rors  influence  the  results  of  the  method.  The  algorithm 
hau  been  formulated  in  ALGOL  and  tested  on  an  IBM  7040  com¬ 
puter;  the  ALGOL  program  is  given. 

Kremer,  H. ,  "Design  of  Linear  Networks  with  Prescribed  Impulse  Rasponse", 
Extended  paper  from  IEEE  Conference  on  Computer-Aided  Design  of  Electronic 
and  Microwave  Circuits  and  Systems,  Hull,  United  Kingdom,  12-14  July  1977. 

A  method  for  the  design  of  linear  networks  with  prescribed 
impulse  response  is  presented,  that  uses  Prony's  algorithm 
for  the  approximation  and  gives  the  network  as  a  parallel 
structure  of  at  most  second  order  subnetworks.  Following 
a  brief  description  of  the  algorithm  some  numerical  and 
computational  aspects  are  discussed.  The  method  is 
illustrated  by  two  examples,  the  design  of  a  Hilbert  filter 
and  an  ideal  low-pass  filter. 

Kulp,  R.C.,  "An  Optimum  Sampling  Procedure  for  Use  with  the  Prony  Method", 
IEEE  Trans,  on  Electomagnetic  Compatibility.  Vol.  EMC-23,  No.  2,  May  1981, 
pp.  67-71. 


The  ill-posed  nature  of  identifying  the  poles  of  a 
transient  response  is  discussed  along  with  a  sampling 
scheme  for  minimizing  the  effect  of  this  ill-posedness . 

Two  simulation  examples  are  given.  The  results  of  the 
second  example  are  compared  to  previously  published 
results. 

Kunz,  K.S.,  et  al.,  "Surface  Current  Injection  Techniques:  A  Theoretical 
Investigation",  IEEE  Transactions  on  Nuclear  Science,  Vol.  NS-2S,  No.  6, 
pp.  1422-1427,  December  1978. 


This  theoretical  study  was  undertaken  to  find  surface 
current  injection  technique  configurations  capable  of  quality 
simulation.  A  configuration  capable  of  quality  simulation 
was  found.  It  compared  favorably  with  conventional  EMP 
simulation  techniques  when  both  were  measured  against 
theoretical  threat  aircraft  responses.  Some  preliminary 
experimental  results  further  confirm  the  capabilities  of 
this  configuration  while  validating  the  theoretical  predic¬ 
tions.  Prony*s  method  was  used  to  find  the  poles  and  residues 
of  the  data. 

Kunz,  K.S.  and  J.F.  Prewitt,  "Practical  Limitations  to  a  Natural  Mode 

Characterization  of  Electromagnetic  Transient  Response  Measurements," 

IEEE  Transactions  on  Antennas  and  Propagation.  Vol.  AP-28,  No.  4,  pp. 

577,  July  1980. 

Practical  limitations  arising  from  limited  slgnal-to-nolse 
ratios  In  actual  translei.t  electromagnetic  measurements 
limits  the  number  of  natural  modes  which  can  be  extracted 
from  such  measurements.  This  ntimber  is  shown  here  to  be 
three  to  six  for  a  complex  scatterer.  In  this  case  an 
aircraft,  five  for  a  "fat"  cylinder  and  thin  wire.  This 
is  fcr  double  exponential  waveform  excitation  of  the  air¬ 
craft  and  step  function  excitation  of  the  fat  cyllrier  and 
thin  wire,  A  few  more  modes  can  be  resolved  if  the  excita¬ 
tion  is  altered  so  as  to  provide  more  high  frequency  energy, 
as  can  be  accomplished  by  changing  the  step  excitation  to  a 
delta  function  excitation,  or  by  signal  conditioning,  for 
example  making  a  derivative  measurement  that  emphasizes  the 
high  frequency  response.  Even  so,  far  fewer  modes  can  be 
resolved  experimentally  than  can  be  predicted.  Experimental 
electromagnetic  transient  response  measurements  should 
therefore  be  made  with  these  practical  limitations  in  mind. 
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Lager,  D.L. ,  H.G.  Hudson,  and  A.J.  Pogglo,  "User's  Manual  for  SEMPEX:  A 
Computer  Code  for  Extracting  Complex  Exponentials  from  a  Time  Waveform" 
AFWL  Mathematics  Note  //AS,  March  1977. 

This  report  Is  the  user's  manual  for  the  SEMPEX  computer 
code,  singularity  Expansion  Method  P^olc  Extraction  Is  a 
technique  In  electromagnetics  In  which  the  free  response 
of  a  structure  Is  expressed  as  a  weighted  sum  of  complex 
exponentials  where  the  damping  factors  of  the  exponentials 
are  derived  from  the  poles  (l.e. ,  singularities)  of  the 
transfer  function  of  the  object. 

Lager,  D.L.,  "Using  Parct-  (Parameter  Estimation  Technique)  for  the  De¬ 
termination  of  Mode  Shapes,"  Lawrence  Livermore  Laboratory,  UCID-17898, 

8  September  1978. 

This  report  applies  the  Parameter  Estimation  Technique 
(PARET)  to  the  measured  responses  of  a  T-shaped  aluminum 
plate.  The  frequency  and  dampln"  for  the  natural  resonances 
and  display  animated  mode  shapes  for  the  plate  are  obtained. 

The  measurements  consist  of  the  response  of  an  accelerometer 
mounted  on  a  comer  of  the  plate  to  hammer  blows  at  38 
separate  locations  on  the  plate.  PARET  is  used  to  compute 
the  natural  resonances  by  determining  the  parameters  of  an 
exponential  model  that  fits  the  responses.  Tabulation  of 
the  frequencies  and  damping  for  seven  modes  occurring  over 
a  frequency  range  SOO-1650  Hz  and  also  plot  the  shape  for 
each  mode  are  made. 

Lager,  O.L.,  "A  Statistical  Evaluation  of  the  Performance  of  the  Time- 
Domain  Paret  Algorithm,"  Lawrence  Livermore  Laboratory,  UClD-18120, 

1  June  1979. 

A  methodology  was  developed  for  evaluating  the  performance 
of  the  time-domain  PARET  (TDP)  algorithm  when  processing 
waveforms  contaminated  with  additive  noise.  The  methodology 
is  general;  it  may  be  applied  to  evaluating  the  performance 
of  other  algorithms  and  provides  a  good  basis  for  comparison. 

The  performance  is  evaluated  by  computing  the  bias,  variability, 
and  root  mean  square  error  in  the  parameters  estimated  by  the 
algorithm  when  processing  ensembles  of  waveforms.  The  wave¬ 
forms  used  in  this  study  were  the  displacements  as  a  function 
of  time  from  a  computer  simulation  (by  the  SAP  IV  Program) 
of  an  eight-floor  structure. 


Lai,  Hang  H.  and  John  D.  Dyson,  "The  Deteimlnatlon  cf  the  Direction  of 
Arrival  cf  an  Interference  Field",  Radio  Science,  Vol.  16,  No.  3, 
pp.  365-376,  May-June  1981. 
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The  problem  of  resolving  the  complex  Interference  field  at 
a  high  receiving  site,  that  is  caused  by  multimode  propaga¬ 
tion  at  high  frequencies,  into  the  individual  plane  wave 
components  and  determining  the  direction  of  arrival  of 
these  components  Is  considered.  The  method  developed  is 
based  on  Prony's  algorithm.  To  Illustrate  the  effect  of 
noise  on  the  solution,  a  first-order  linearized  model  of  a 
statistical  study  of  a  two  plane  wave  problem  for  this 
conventional  Prony's  method  Is  presented.  A  modification 
of  Prony's  technique,  an  eigenvalue  approach,  allows  one 
to  analytically  determine  the  number  of  waves  in  the 
problem,  and  this  approach  is  lass  susceptible  to  the 
effect  of  noise.  Data  obtained  from  a  controlled  experi¬ 
ment  when  a  signal  from  a  distant  source  was  repeated  and 
retransmitted  at  a  different  location  verify  that  the 
system  as  wall  as  the  algorithm  Is  capable  of  resolving 
waves  coming  from  different  directions.  Experimental 
evidence  related  to  multimode  as  well  as  single-mode 
propagation  coming  from  the  same  source  has  been  obtained. 

Leveckls,  A.S.,  "Application  of  Singularity  Expansion  to  Electromagnetic 

Pulse  Response  Waveforms,"  Master's  Thesis,  Massachusetts  Institute  of 

Technology,  June  1977. 

Here  la  presented  the  results  of  an  attempt  to  use  Prony's 
method  on  EMF  response  waveforms  on  ships.  The  technique 
to  acquire  the  experimental  data  la  detailed;  the  numerical 
problem  Is  posed,  and  the  goals  are  quickly  simplified  to 
attempting  to  find  the  best  curve-fit  to  the  waveform,  using 
complex  exponentials,  and  then  hoping  to  eliminate  the  curve- 
fit  poles.  A  method  Is  presented  which  aids  In  the  'physical 
pole'  selection  process,  and  explained  also  Is  the  signi¬ 
ficance  of  the  right  half  plane  poles  In  the  results. 

Lytle,  R.J.  and  D.L.  Lager,  "Using  the  Natural-Frequency  Concept  In  Re¬ 
mote  Probing  of  the  Earth,"  Radio  Science.  Vol.  11,  No.  3,  pp.  199-209, 

March  1976. 

A  technique  Is  given  for  determining  the  thicknesses  and 
the  electrical  constitutive  parameters  of  a  planar-layered 
medium  such  as  .i  coal  seam  In  a  mine  environment.  Time- 
domain  experimental  data  are  analyzed  with  Prony's  method 
to  determine  the  natural  frequencies  of  the  layered  medium. 
Explicit  relations  are  given  (for  dielectric  layers)  for 
determining  the  thicknesses  and  dielectric  constants  from 
the  experimentally  determined  natural  frequency  results. 

Explicit  expressions  are  also  given  (for  conductive  layers) 
for  calculating  the  electrical  thickness  from  natural  fre¬ 
quency  results.  The  method  Is  Illustrated  with  sample 
numerical  results. 
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Makhoul,  J. ,  "Linear  Prediction:  A  Tutorial  Review",  Proceedings  of  the 
IEEE.  Vol.  63,  No.  4,  pp.  561-580,  April  1975. 

This  paper  gives  an  exposition  of  linear  prediction  in  the 
analysis  of  discrete  signals.  The  signal  is  modeled  as  a 
linear  combination  of  Its  past  values  and  present  and  past 
values  of  a  hypothetical  input  to  a  syater  whose  output  is 
the  given  signal.  In  the  frequency  domain,  this  is  equiva¬ 
lent  to  modeling  the  signal  spectrum  by  a  pole-zero  spectrum. 

The  major  part  of  the  paper  is  devoted  to  all-pole  models. 

The  model  parameters  are  obtained  by  a  least  squares  analysis 
in  the  time  domain.  Two  methods  result,  depending  on  whether 
the  signal  is  assumed  to  be  stationary  or  nonstationary.  The 
same  results  are  then  derived  In  the  frequency  domain.  The 
resulting  spectral  matching  formulation  allows  for  the  modeling 
of  selected  portions  of  a  spectrum,  for  arbitrary  spectral 
shaping  In  the  frequency  domain,  and  for  the  modeling  of 
continuous  as  well  as  discrete  spectra.  This  also  leads  to 
a  discussion  of  the  advantages  and  disadvantages  of  the 
least  squares  error  criterion.  A  spectral  interpretation  is 
given  to  the  normalized  minimum  prediction  error.  Applica¬ 
tions  of  the  normalized  error  are  given,  including  the  de¬ 
termination  of  an  "optimal"  number  of  poles.  The  use  of 
linear  prediction  in  data  compression  is  reviewed.  For 
purposes  of  transmission,  particular  attention  is  given  to 
the  quantization  and  encoding  of  the  reflection  (or  partial 
correlation)  coefficients.  Finally,  a  brief  introduction 
to  pole-zero  modeling  is  given. 

i-iarla,  G.A.  and  M.M.  Fahmy,  "Ip  Approximation  by  Exponentials",  IEEE 
Tranaactlons  on  Circuit  Theory,  pp.  71-74,  January  1973. 

An  interatlve  method  Is  proposed  for  obtaining  a  solution 
to  the  problem  of  1*°  approximation  by  exponentials.  The 
method  always  converges  and  has  the  advantage  that  the 
samples  need  not  be  equidistant.  Using  the  proposed  method 
and  allowing  p  to  Increase  yields,  in  a  much  easier  way, 
results  that  are  comparable  to  those  obtained  by  Flschl  for 
the  Chebyshev  approximation. 

Marzollo,  A.,  "On  the  Mean  Square  Approximation  of  a  Function  with  a 
Linear  Combination  of  Exponentials",  Int.  Journal  Control.  Vol.  9,  No.  1, 
pp,  17-26,  1969. 

This  paper  deals  with  the  problem  of  finding  the  2N 
coefficients  of  a  linear  combination  of  N  exponentials 
which  minimizes  Its  mean  square  error  from  a  given 
function  f(t).  The  real  difficulty  lies  In  the  determina¬ 
tion  of  the  optimal  exponents. 
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A  raal  function  of  N  rtal  varlablas  li  conatructad,  auch  chat 
cha  polnta  whara  It  ia  minimal  giva  tha  optimal  axponanta. 

Such  a  function  may  ba  conaidarad  aa  tha  aampla  for  t  T  of 
cha  raaponaa  to  f(T  c)  of  a  aimpla  ayacam.  Thia  implamanta- 
cion  allowa  a  aignificant  aimplif icacion  of  cha  givan  problam 
from  a  computational  point  of  viaw,  for  inatanca  by  hybrid 
machoda. 

Millar,  E.K. ,  at  al.,  "Radar  Targat  Claaalficatlon  Ualng  Tamporal-Moda 
Analyaia",  Lawranca  Livarmora  Laboratory,  UCRL-5182S,  27  May  1975. 

Exploratory  taat  raaulta  ara  pratantad  for  a  naw  radar 
targat-clasaificrition  tachniqua  baaad  on  procaaalng  tha 
transiant  acattorad  wavaform  from  a  targat  uaJng  tamporal- 
ffloda  analyaiii.  Tha  naw  mathod  offara  tha  potaiTitial 
advantagaa  of  targat  charaetarlzatlon  indapondant  of 
aapact  anglt  and  polarization,  a  maana  for  quicldy 
comparing  targata  againat  a  atorad  library,  and 
aubatantially  aimplif ying  data  storaga  and  procaaaing 
raquiramanta.  Although  aavaral  thaoratical  quaationa 
ramain  to  ba  anawarad,  initial  taat  raaulta  look  quita 
promiaing  for  thia  naw  approach. 

Millar,  E.K. ,  J.N.  Brittingham,  and  J.L.  Willowa,  "Tha  Darivation  of 
Simpla  Polaa  in  a  Tranafar  Function  from  Raal  Fraquancy  Information 
Fart  3:  Ob j act  Claaaification  and  Idantif ication" ,  Lavrranca  Livarmora 
Laboratory,  UCRL»S2211,  10  January  1977. 

Thia  raport  ia  tha  third  in  a  aariaa  of  thraa  that 
avaluata  a  tachniqua  (fraquancy-domain  Prony)  for 
obtaining  tha  polaa  of  a  tranafar  function.  Tha  main 
objactiva  waa  to  aaaaaa  tha  faaaibility  of  claaaifying  or 
idantifying  ahip-llka  targata  by  uaing  pola  aata  darivad 
from  fraquancy-domain  data.  A  pradictor-corralator 
procadura  for  uaing  apactral  data  and  library  pola  aata 
for  thia  purpoaa  waa  davalopad.  Alao  atudiad  waa  an 
itarativa  mathod  for  raducing  tha  af facta  of  noiaa  and  a 
tachniqua  baaad  upon  magnituda-only  apactral  data. 

Miller,  E.K. ,  "Data  Charaetarlzatlon  and  Compraaaion" ,  Lawranca  Livarmora 
Laboratory,  UCID-17511,  6  July  1977. 

Tha  information- tranaformation  procaaa  baaad  upon  Prony 'a 
mathod  haa  baan  found  to  ba  Incraaalngly  uaaful  for  appli¬ 
cation  to  elactromagnatlc-typa  problama  in  particular  and 
a  growing  variety  of  phyaical  problama  in  ganaral. 

Millar,  E.K.  and  D.L.  Lagar,  "Information  Extraction  Uaing  Prony'a  Mathod" 
Lawranca  Livarmora  Laboratory,  UCRL-52329,  23  Auguat  1977. 


This  raport  consldara  tha  quaatlona  of  information  contant 
and  data  tranaformatlon.  To  daflna  Information  contant, 
tha  concapta  of  data  rank  (numbar  of  polaa  in  tha  axponantial 
data)  and  data  pracialon  ara  uaad.  To  atudy  data  trana- 
formation,  a  vavaform  la  dafinad  in  karma  of  apacifiad 
polaa  and  raalduaa,  ia  aamplad  to  obtain  tha  aaaoclatad 
data,  and  ia  oparatad  on  ualng  Prony'a  mathod  to  tranaform 
tha  data  and  compara  it  with  tha  original  information 
contant.  Othar  numarical  procaaaaa,  auch  aa  matrix 
invaralon,  can  ba  vlawad  in  tha  aama  way*  Tha  problama  of 
ill-conditioning  ara  ralatad  mora  to  tha  data  than  to  tha 
procaaa.  Finally,  applicationa  of  Prony'a  mathod  to  othar 
phyaical  problama  with  axponantial  aolutiona  are  conaldarad. 

Miller,  E.K.  and  O.L.  Lagar,  "Prony'a  Mathod  for  tha  Angla  Domain", 
Lawranca  Livarmora  Laboratory,  UC1D-17S02,  Rav.  1,  6  October  1977. 

Prony'a  mathod,  a  tachniqua  for  obkainlng  tha  paramatara 
of  a  aariaa  of  complax  axponantiala ,  ia  uaad  hara  to 
analyze  tha  radiation  pattern  of  a  linear  array  of  lao- 
tropic  radiatora.  Tha  locationa  and  amplitudaa  of  the 
radlatora  ara  derivable  from  tha  pattern.  Poaaibla 
applicationa  include  imaging  and  array  aynthaaia. 

Millar,  E.K.  and  D.L.  Lagar,  "Radiation  Field  Analyaia  and  Synthaaia  Uaing 
Prony'a  Mathod",  Lawranca  Livarmora  Laboratory,  UCRL-80532,  16  January 
1978. 


Prony'a  mathod,  a  tachniqua  for  extracting  tha  paramatara 
of  an  axponantial  aariaa  from  Ita  diacrata  aamplaa,  ia 
applied  hara  to  far-fiald  pattern  data.  Tha  mathod  offara 
poaalbilltiaa  in  imaging  aourca  diatributlona  and  in  array 
analyaia  and  aynthaaia. 

Millar,  E.K. ,  "Prony'a  Mathod  Ravlaitad",  Lawranca  Livarmora  Laboratory, 

UCRL-52590,  18  October  19*’'. 

Thia  report  continuaa  evaluation  of  tha  moving-win  ^ 
approach  for  reducing  tha  dalatarioua  affacta  if  •■>..ae  to 
which  Prony  procaaaing  ia  particularly  aanaitlva.  a 
aariaa  of  aystamatic  computer  axparlmanta  have  bean 
conducted  on  analytically  apacifiad  tranaiant  wavaforma 
and,  in  particular,  have  axaminad  tranamiaaion  of  arrora 
though  tha  pmcaaa  by  determining  tha  accuracy  of  tha 
input  data,  tha  charactariatic-aquation  coefficlanta  and 
roota,  and  tha  polaa.  Statlatical  diatributlona  of  thaaa 
quantitlaa  for  comblnatlona  of  input  data  paramatara  have 
bean  derived. 


E.K.  and  D.L.  Lagar,  "Invaralon  of  Ona-Dlmanaional  Scattaring 
Data  Ualng  Prony'a  Mathod",  Lavranca  Livaraora  Laboratory,  UCRL-52667, 

12  Fabruary  1979. 

A  ona-dlaanalonal  configuration  it  tha  almplaat  gaonatry 
to  invar t,  yat  it  haa  practical  application  to  auch 
problama  aa  acattaring  from  inhonoganao m  half-apacaa  and 
propagation  on  non-uniform  tranamiaalon  linaa.  Whathar 
tha  madiua  paramatara  vary  contlnuoualy  or  diacrataly  with 
poaltion,  tha  problam'a  numarical  daacription  can  uaually 
ba  davalopad  in  finlta-tima  approximation.  Aa  auch,  tha 
acattarad  and  tranamittad  flalda  can  ba  rapraaantad  aa 
axponantlal  aariaa,  whoaa  axponanta  ara  ralatad  to  tha 
alactrical  thicknaaaaa  of  tha  lay ara  which  maka  up  tha 
modal.  If  tha  axponanta  or  polaa  ara  darivabla  from  fiald 
data,  than  tha  invaraa  problam  ia  formally  aolvabla.  Thla 
raport  conaidara  application  of  Prony'a  mathod,  a  procadura  , 
for  obtaining  tha  polaa  of  axponantlal  aignala,  to  auch 
ona-dimanaional  problama.  Analyaia  of  both  tima-domain  and 
fraquancy-domain  data  ia  atudiad.  Tha  af facta  of  tha 
madium  charactariatica,  numbar  of  layara,  and  othar  factora 
ara  axaminad.  It  ia  concluded  that  Prony'a  mathod  haa 
marit  for  cartain  claaaaa  of  ona-dimanaional  invaraa 
problama. 

Millar,  E.K. •  "A  Study  of  Targat  Zdantiflcation  Uaing  Blactromagnatic 
Polaa",  Lawranca  Llvarmora  Laboratory,  UC11L-S268S,  IS  March  1979. 

Thraa  tachnlquaa  to  identify  radar  targata  by  ualng  thalr 
alactromagnatic  raaonaneaa  or  polaa  ara  atudiad.  A 
aingular  advantage  of  auch  a  pola-baaad  approach  ia  that 
it  can  operate  indapandantly  of  tha  targat  excitation.  A 
library  of  pole  aata  ia  uaad  both  to  generate  tha  targat 
wavaforma  (to  which  nolaa  ia  added)  and  to  correlate  with 
paramatara  extracted  from  thoaa  wavaforma.  Tha  tachnlquaa 
atudiad  ara  linear  prediction,  raaldua  calculation,  and 
pole  calculation.  Direct  waveform  correlation  waa  alao 
evaluated  for  compariaon  with  tha  pola-baaad  tachnlquaa. 

Uaing  10- run- average  correlation  valuaa,  arror-froa 
identification  waa  accompllahad  by  pole  calculation  at  a 
paak-aignal-to-paak-noiaa  ratio  aa  low  aa  10,  and  by  wave¬ 
form  correlation  at  a  nolaa  ratio  of  unity. 

Millar,  E.K. ,  D.L.  Lagar  and  J.t.  Okada,  "Imaging  of  Linear  Source 
Diatributlona" ,  Lawranca  Livermore  Laboratory,  UCRL-52822,  Movambar  1979. 

An  approach  baaad  on  Prony'a  mathod  ia  uaad  to  find  tha 
apatial  polaa  of  atraight  wlraa  of  varying  length  and 
radiua  from  the  far-fiald  pattarna.  For  the  antenna  caaa, 
polaa  having  tha  largaat  raalduaa  ara  found  near  tha  teed 
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region  and  antanna  anda;  amallar  onaa  ara  found  outalda 
chaaa  raglonii  Tha  polau  ara  attributed  to  radiation  from 
the  exciting  aourca  and  wire  anda,  due  to  charge  accelera¬ 
tion,  and  from  a  ailghtly  damped  traveling  wave.  In  tha 
acattarlng  caaa,  the  polaa  with  tha  largaat  reaiduaa  ara 
located  at  tha  wire's  anda  for  naar-axlal  Incidence  and 
ara  again  dua  to  and  radiation.  By  examining  tha  residua 
magnitudes  as  a  function  of  wire  radius,  a  correlation  is 
found  with  radiation  from  tha  feed  region  of  an  Infinite 
antanna  and  from  tha  anda  of  a  saml-lnflnlta  antanna. 

Millar,  G.,  "Representation  and  Analysis  of  Signals;  Part  XXVI,  Least- 
Squares  Approximation  of  Functions  by  Exponentials",  Dspartmant  of 
Electrical  Engineering,  Tha  John  Hopkins  University,  Baltimore,  Maryland, 
June  1969. 

The  approximation  of  an  analytic  time  function  In  the  least- 
squares  sense  by  sums  of  exponentials  Is  considered  from 
several  different  points  of  view.  McDonough's  method  for 
finding  the  poles  Is  derived  by  three  different  approaches. 

In  the  process,  a  new  method  la  developed  which  offers  the 
advantages  of  the  earlier  resulto  achieved  by  McDonough  and 
by  McBride,  Schaefgen,  and  Stelglltx.  This  new  method 
reveals  the  link  between  chase  earlier  methods  and  provides 
a  standard  for  comparing  chase  two  linear  iterative  schemas 
using  several  numerical  examples. 

Mlctra,  R.  and  C.O.  Stearns,  "Field-Strength  Measurements  In  a  Multipath 
Field  Using  Linear  and  Circular  Proving",  Radio  Science.  Vol.  2,  Mo.  1, 
pp.  101-110,  January  1967. 

Techniques  for  resolving  Che  componeuts  of  a  mu.lcipach  field 
are  investigated  In  this  paper.  Both  linear  and  circular 
probing  schemes  are  discussed.  A  method  for  analysing  Che 
linearly  probed  data  Is  given,  which  Is  consldtrably  simpler 
Chan  the  procedure  followed  by  Watterson  In  an  «Hi.ller  paper. 

The  method  la  alio  extended  to  the  circular  case.  It  Is 
found  Chat  the  ttisoluclon  In  the  circular  case  Is  poorer  Chan 
for  Che  linear  scheme. 

Moffett,  D.L.  and  K.A.  Shubert,  "Natural  Reionancas  Via  Rational 
Approximations",  IEEE  Transactions  on  Antennas  ar.g  Propagation.  Vol.  AP-25, 
No.  3,  pp.  657-660,  September  1977. 

An  approach  for  locating  the  dominant  complex  natural 
resonances  (poles)  of  scatt^rers  or  radiators  Is  described. 

Admittance  valuee,  calibrated  using  Richmond's  thin  wire 
comouter  programs,  at  a  regular  grid  of  complex  and  real 
frequencies  are  fit  by  a  rational  function.  Zeros  of  the 
resultant  numerator  and  denominator  polynomials  then  yield. 


raipactivaly,  the  ztroi  and  polaa  of  tho  atructura.  Tha 
mathod  la  llluatratad  for  raaonancaa  of  atralghc  and  bant 
wlra  atructuraa. 

Moffatt»  D.L.,  L.C.  Chan  and  G.A.  Huwiahar,  "Charactarizatlon  of 
Subaurfaca  Elactromagnatlc  Soundlnga"»  Ohio  Stata  SlactroSclanca 
Laboratory,  Annual  Tachnlcal  Raport  4A90-1,  Saptambar  1977. 

In  thla  raport,  datalla  on  a  flrat  yaar'a  afforta  to 
charactariaa  aubaurfaca  radar  targat  algnaturaa  arw 
aummariaad.  Tha  raport  covara  thraa  main  toplca;  an 
Initial  maeauramant  program  to  obtain  tranalant  radar 
algnaturaa  of  burlad  objacta,  a  daacriptlon  of  a  compra- 
hansiva  computar  program  to  axtract,  by  varioua  maana,  tha 
complax  natural  raaonancaa  of  burlad  objacta,  and  tha 
complicatlona  involved  in  axtanding  a  aynthatic  radar 
imagery  technique  to  burlad  targata.  In  a  concluaiona 
aaction,  tha  atata  of  tha  art  la  aummariaad  and  raaaarch 
for  tha  next  yaar'a  effort  ia  diacuaaad. 

Moffatt,  D.L.  and  L.C.  Chan,  "Charactariaatlon  of  Subaurfaca  Soundinga", 
Ohio  Stata  ElactroScianca  Laboratory,  Final  Raport  4490-2,  January  1979. 

Raaaarch  on  a  National  Science  Foundation  Grant  for  tha 
period  1  Juna  1976  to  31  Dacambar  1978  ia  aummariaad. 

Unique  new  axcltation  invariant  paramatara  (complex 
natural  raaonancaa)  of  aubaurfaca  anomallaa  are  daacribad. 

Mathoda  for  extracting  and  exploiting  thaaa  excitation 
invariant  paramatara  in  tha  praaanca  of  noiaa  and  clutter 
are  given.  Stata-of-tha-art  aubaurfaca  radar  ayatama  are 
illuatratad.  Uaing  auch  a  ayatam,  target  charaetariaation 
procaduraa  era  auccaaafully  illuatratad  for  a  variety  of 
subaurfaca  targata  using  real  radar  data. 

Moffatt,  O.L.  ft  al,  "Transient  Raaponaa  Charactariatics  in  Idantificctlon 
and  Imaging",  IEEE  Tranaactiona  on  Antannaa  and  Propaaation.  Vol.  AP-29 , 
No.  2,  pp.  192-20S,  March  1981. 

Thraa  differuut  soluLiona  to  problems  in  identification  and 
invatsa  scattaring  using  raatrictad  far-xona  acattarlng 
data  are  ravlavad  and  illuatratad.  The  data  are  raatrictad 
in  tha  sense  that  tha  fraquanclao,  aspects,  and  polarlxations 
ara  liraitad.  All  throe  solutions  are  basad  in  part  on  a 
fundamental  time  domain  viewpoint  wharaby  tha  far-rona 
scattering  characteilr tica  of  any  finita  object  ara  uniquely 
aumnarixad  by  its  impulse  raaponaa  waveforms.  Insight 
provided  by  these  wavatorma  ia  exploited  via  geometrical 
charactariatics  extracted  from  ataady-atata,  forced,  <*nd 
free  rasponsaa  of  tha  object.  For  tha  present  purposes  tha 
fundamental  importance  of  signaling  waveforms  whose 


wtvalttngths  art  within  an  ordar  of  magnltuda  of  tha  objact 
dlmanalon  la  damonatratad.  Tha  baalc  mathoda  dlacuaaad 
hava  baan  pravloualy  Idantlflad  in  tha  lltaratura  aa 
natural  raaonanca  aatlmation,  radar  Imaging  from  ramp 
raaponaa  algnaturaa,  and  low  fraquancy  claaaiflcatlon. 

Ralatad  common  fatturaa,  additional  Inalght,  and  aoma  naw 
raaulta  ara  glvan.  Applicatlona  of  all  thraa  mathoda  to 
objacta  ranging  from  almpla  to  complax  gaomatrlaa  ara 
daacrlbad. 

Murphy,  W.D.  and  M.L.  Bamaba,  "Mumarlcal  Procaduraa  for  Solving 

Monaymmatrlc  Eiganvalua  Problama  Aaaoclatad  with  Optical  Raaonatora'* 

Appliad  Optica.  Voli  17,  No.  15,  pp.  2358-2365,  1  Auguat  1978. 

Tha  Prony  mathod  la  axtandad  to  handla  tha  nonayrmnatrlc 
algabralc  alganvalua  problam  and  improvad  to  aaarch  auto¬ 
matically  for  tha  numbar  of  dominant  alganvaluaa.  A 
alngla  Itaratlva  algorithm  la  glvan  to  computa  tha 
aaaoclatad  alganvactora.  Raaolutlon  atudlaa  ualng  tha  QR 
mathod  ara  mada  In  ordar  to  datarmlna  tha  accuracy  of  tha 
matrlc  approximation.  Numarlcal  raaulta  ara  glvan  for  both 
almpla  wall-daflnad  raaonatora  and  mora  complax  advancad 
daalgna  containing  multlpla  propagation  gaomatrlaa  and 
mlaallgnad  mlrrora. 


McDonough,  R.N. ,  "Rtpras«nt«tlon  and  Analysis  of  Signals,  Part  XV:  Matchad 
Exponents  for  ths  Rsprssantatlon  of  Signals,"  Ths  Johns  Hopkins  University, 
Department  of  Electrical  Engineering,  Baltimore  Maryland,  30  April  1963. 

Approximation  by  sums  of  exponentials  is  considered  from 
various  points  of  view,  such  as  the  time-domain  approximation 
problem  of  network  synthesis,  In  which  an  analytic  time 
function  Is  to  be  approximated,  or  situations  In  which  dis¬ 
crete  measured  values  of  some  dynamic  process  are  to  be 
smoothed  and  fitted  using  an  exponential  model.  Prony's 
method  for  approximate  determination  of  the  exponents  Is 
exhaustively  reviewed  In  all  Its  forms. 

McDonough,  R.N.  and  W.H.  Huggins,  "Best  Least-Squares  Repressntatlon  of 
Signals  by  Exponentials,"  IEEE  Transactions  on  Automatic  Control.  Vol. 
AC-13,  No.  4,  pp.  408-412,  August  1968. 

In  this  paper  the  approximation  of  a  given  real  time  function 
over  (0,«*)  by  a  linear  combination  of  a  glvan  number  n  of 
exponentials  Is  considered,  such  that  the  Integrated  squared 
error  la  mlnlmlzad  over  both  the  n  coefficients  of  the  linear 
combination  and  the  n  exponents  used.  The  usual  nacassary 
condition  for  statlonarlty  of  the  Integrated  squared  error 
leads  to  a  set  of  2n  simultaneous  equations,  nonlinear  In  the 
exponents.  This  condition  la  Interpretsd  In  the  geometric 
language  of  abstract  vsetor  spaces,  and  an  equivalent  con¬ 
dition  Involving  only  the  exponents,  with  the  coefficients 
suppressed,  Is  dsveloped.  It  Is  next  Indicated  how  this 
latter  condition  can  be  applied  to  signals  which  ara  not 
known  analytically,  but  only,  for  exampls,  as  voltages  re¬ 
corded  on  magnetic  tape,  or  as  a  tabls  of  samplsd  values. 

The  condition  still  In  effect  requires  solution  of  nonlinear 
algebraic  equations,  and  a  linear  Iterative  method  Is  pro¬ 
posed  for  this  purpose.  Finally,  the  procedure  Is  Illustrated 
with  a  simple  example. 


Osborne,  M.R.,  "A  Class  of  Nonlinear  Regression  Problems",  in  Data 
Representation.  University  of  Queensland  Press,  pp.  94-100,  1970. 


This  paper  considers  the  problem  of  determining  the  para¬ 
meters  a,  8  in  the  model  equation 

M*  •• 
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y(t)  "2  ®i  ^k  **•» 


from  observation 
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made  at  a  sequence  of  distinct  time  point  t. •  t.  ...,  t  • 
It  is  assumed  that  the  e.  (the  experlmental^errofs)  are 
independent,  normally  distributed,  random  variables  with 
mean  zero  and  standard  deviation  a  so  that  least  squares 
methods  are  appropriate.  It  is  also  convenient  to  assume 
that  the  functions  f .(t,  8)  form  a  Chebyshev  set  for  the 
values  of  8  of  Interest.  “  This  means  that  no  nontrivial 
linear  comSination  of  the  f  can  vanish  identically  on  any 
set  of  k  distinct  values  of^t.  It  is  convenient  to  refer 
to  a  regression  problem  of  the  above  form  as  separable,  and 
refer  to  a  and  8  as  the  linear  and  nonlinear  parameters 
respectively. 
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Parsons,  D.H. ,  "Biological  Problems  Involving  Sums  of  Exponential 
Functions  of  Time:  A  Mathematical  Analysis  that  Reduces  Exparlmental  Time" 
Mathematical  Blosclences .  #2,  pp.  123-128,  1968. 

The  analysis  given  In  this  article,  which  Is  based  on  the 
theory  of  difference  equations,  enables  one  to  determine 
the  Integer  rank  In  advance,  and  then  determine  the  2p 
unknown  constants  by  means  of  simultaneous  linear  equations. 

Paul,  A.K. ,  "Anharmonlc  Frequency  Analysis",  Mathematics  of  Computation. 
Vol.  26,  No.  118,  pp.  437-447,  April  1972. 

A  new  numerical  method  of  frequency  analysis  Is  described, 
designed  mainly  to  search  for  discrete  frequencies  In  a 
time  series.  An  Integral  transform  Is  applied  twice  to  the 
data  for  different  reference  tines.  A  complex  amplitude 
within  a  selected  narrow  frequency  band  is  obtained  for  each 
transform.  The  frequency  la  then  determined  from  the  phase 
change  of  the  complex  amplitude  over  the  different  of  the 
two  reference  times.  Very  high  precision  Is  obtained,  which 
Is  demonstrated  In  two  examples. 

Pearson,  L.W. ,  M.L.  Van  Blarlcum,  and  R.  Mlttra,  "A  New  Method  for  Radar 
Target  Recognition  Based  on  the  Singularity  Expansion  for  the  Target", 
Proceedings  of  the  1975  IEEE  International  Radar  Conference.  Arlington, 
Virginia,  pp.  452-456,  21-23  April  1975. 

The  results  of  an  Initial  study  of  the  use  of  the  recently 
Introduced  Singularity  Expansion  Method  (SEM)  characterization 
of  a  scatterer  as  a  means  of  radar  target  Identification  are 
reported.  Three  significant  advantages  of  this  method  over 
target  Identification  schemes  are  explored.  These  are  1)  the 
aspect-dependence  of  transient  backscattered  waveforms  can  be 
suppressed,  2)  the  transmitted  waveform  requirements  are  only 
moderately  demanding  on  hardware,  and  3)  the  SEM  poles  consti¬ 
tute  a  smaller  collection  of  numbers  to  be  presented  to 
pattern  recognition  algorithms.  Examples  are  presented  Indi¬ 
cating  the  following  features  of  the  method:  viability  In  the 
presence  of  noise;  excitability  of  a  reasonable  number  of 
natural  frequencies  with  a  moderate  bandwidth  pulse;  and  the 
ability  of  a  few  poles  to  Identify  a  target  from  among  a  col¬ 
lection  of  known  targets. 

Pearson,  L.W.  and  D.  R.  Roberson,  "The  Extraction  of  the  Singularity 
Expansion  Description  of  a  Scatterer  from  Sampled  Transient  Surface 
Current  Response",  IEEE  Transactions  on  Antennae  and  Propagation. 

Vol.  AP-28,  No.  2,  pp.  182-190,  March  1980. 


3-32 


■V  m  ' 


A  method  Is  presented  whereby  one  cen  extract  the  singu¬ 
larity  expansion  method  (SEM)  description  o£  an  object's 
electromagnetic  scattering  response  from  spatially  sampled 
transient  surface  currents.  The  currents  are  excited  by  a 
known  excitation.  The  SEM  data  are  recoverable  to  the 
degree  chat  spatial  coupling  and  (frequency)  spectral 
Intensity  excite  a  given  SEM  mode.  Results  of  a  numerical 
study  of  the  method  using  Che  transient  response  of  a  Chin 
wire  are  reported.  The  data  used  In  the  study  were 
computed  using  a  time  domain  Integral  equation  technique. 

The  ultimate  utility  of  the  method  lies  In  the  recovery  of 
SEM  data  from  measured  data,  thus  admitting  complex-shaped 
objects  Into  the  realm  of  SEM  description.  The  method  Is 
based  on  e  Prony-type  pole/resldue  extraction  procedure. 

Pearson,  L.W.  and  O.R.  Roberson,  "Extrapolation  of  Induced  Transient 
Surface  Currents  on  a  Scatterer  to  General  Band-Limited  Excitation 
Conditions",  University  of  Kentucky  Electromagnetics  Research  Report, 
pp.  1-49,  October  1978. 

A  technlqixe  for  the  extraction  of  an  approximate  singularity 
expansion  description  of  the  electromagnetic  response  of  a 
scatterer  la  presented.  The  description  Is  extracted  from 
transient  waveforms  which  are  spatial  samples  of  the 
transient  current  on  the  scatterer  excited  by  a  known 
excitation.  The  method  can  use  any  reliable  Prony-type 
processor  for  exponential  analysis  of  the  transient  wave¬ 
forms  as  long  as  It  provides  physlcally-meanlngful  polos 
from  the  waveforms.  The  fact  that  singularity  expansion 
poles  are  common  to  all  waveforms  on  the  object  provides  a 
computational  advantage  over  a  waveform-by-waveform  Prony- 
type  analysis.  A  "consensus  pole  set"  procedure  Is 
described  to  exploit  this  redundancy.  Some  Issues  relating 
to  the  completeness  of  a  set  of  transient  data  to  the  class 
of  excitations  to  which  the  singularity  expansion  Is  applied 
Is  discussed.  The  one-dlmenslonel  example  of  a  thln-wlre 
scatterer  Is  used  to  test  the  practicability  of  the  method. 

The  ultimate  utility  of  the  method  Is  demonstrated  by  way 
of  expansion  of  the  extracted  data  to  new  excitations. 
Recommendations  are  made  toward  an  Implementation  of  the 
method  using  measured  data. 

Perkins,  J.F.  and  Charles  Cason,  "Effects  of  Small  Misalignments  In  Empty 
Unstable  Resonators",  Applied  Physics  letters.  Vol.  31,  No.  3,  pp.  198- 
200,  1  August  1977. 

Effects  of  misalignment  of  rectangular-cross-sectlon 
unstable  resonators  are  Investigated  by  calculations  of 
eigenvalues  and  shapes  of  transverse  modes  of  empty  strip 
resonators.  Small  misalignments  can  substantially  affect 


mode-loss  separation.  A  design  criterion  Is  suggested  for 
.choice  of  Fresnel  numbers  such  that  effects  of  misalign¬ 
ments  are  mlnlm1l,2ed . 

Perl,  W. ,  "A  Method  for  Curve-Fitting  by  Exponential  Functions", 
International  Journal  of  Applied  Radiation  and  Isotropes .  Vol.  8,  pp.  221- 
222. 

A  method  la  presented  for  fitting  exponential  functions  to 
data  of  the  decay  type.  In  an  elementary  "peellng-off" 
form  of  the  method,  the  local  ordinate  Is  plotted  against 
the  local  first  derivative  of  the  data.  The  straight  line 
portion  of  this  graph  at  large  time,  If  It  exists,  yields 
simultaneously  a  constant  (asymptotic  value)  plus  an 
exponential,  which  can  be  subtracted  from  the  given  cur.'e 
and  the  procedvire  repeated.  The  method  permits  judging  the 
data  In  simple  fashion  for  suitability  of  fit  by 
exponentials  simultaneously  with  the  fitting  Itself.  The 
method  Is  formulated  as  an  analytic  version  of  Prony's 
method,  which  permits  clarification  of  certain  difficulties 
of  the  latter.  Procedures  are  suggested  for  fitting 
exponentials  with  closely  spaced  half-lives,  and  for 
fitting  functions  other  than  exponential.  The  method  Is 
Illustrated  by  application  to  data  on  the  time  course  of 
respiratory  uptake  of  a  highly  fat-soluble  Inert  gas  by 
human  subjects. 

Pogglo,  A.J.,  M.L.  Van  Blarlcum,  E.K.  Miller,  and  R.  Mlttra,  "Evaluation 
of  a  Processing  Technique  for  Transient  Data'',  IEEE  Transactions  on 
Antennas  and  Propagation.  Vol.  AP-26,  No.  1,  pp.  165-173,  January  1978. 

A  data  processing  technique,  namely  Prony's  method.  Is 
evaluated.  The  use  of  the  algorithm  for  EMP  and  other 
transient  problems  is  illustrated  and  the  relationship 
between  the  waveform  parameterization  and  the  singularity 
expansion  method  (SEM)  is  described.  The  impulse  response 
of  a  synthesizer  network  is  obtained  using  Prony  processing 
of  the  output  of  the  double  exponential  excited  circuit. 

Difficulties  which  arise  in  this  type  of  processing  such 
as  rank  deficiency,  aliasing,  and  noise  effects  are  con¬ 
sidered  and  methods  for  alleviation  such  as  fllterlug  and 
rank-over specif ication  are  introduced  and  evaluated.  The 
overall  acatus  of  waveform  parameterization  as  applied  to 
SEM  is  considered  and  suggestions  for  future  research  pre¬ 
sented  . 

Prewitt,  J.F.  and  K. S.  Kunz,  "SEM  Pole  Extraction:  A  Feasibility  Study", 
Mission  Research  Corp.,  AMRC-N-&1,  March  1978. 
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This  communication  establishes  typical  experimental 
electromagnetic  transient  response  measurement  limitations, 
and  applies  these  limitations  to  the  natural  mode 
decompoaltlon  of  the  response  of  a  variety  of  scatterers 
due  to  EMP. 

Price,  H.J.,  "An  Improved  Prony  Algorithm  for  Exponential  Analysis",  AFWL 
Mathematics  Notes  #59,  November  1978. 

An  Improved  Prony  algorithm  is  described  which  produces 
a  lower  sum  of  squares  of  residuals  than  the  usual  Prony 
algorithm.  Indeed,  the  new  algorithm  produces  the  lowest 
sun  of  squares  of  residuals  possible.  A  simple  numerical 
example  Is  used  to  compare  the  new  and  old  algorithms. 

Within  the  context  of  the  new  algorithm,  a  technique  to 
constrain  poles  Is  also  described. 

Price,  M.L.,  et  al.,  "Electrical  and  Photon  Tests  of  a  Resonant  Satellite 
Shape,"  IEEE  Transactions  on  Nuclear  Science.  Vol.  NS-25,  No.  6,  pp.  1358- 
1364,  December  1978. 

A  resonant  body  has  been  tested  and  analyzed  to  determine 
Its  SGEMF  response.  Both  electrical  and  photon  excitation 
were  used  and  their  results  compare  favorably  with  pre¬ 
dictions.  A  modal  analysis  of  the  experimental  data  using 
Prony's  algorithm  yielded  the  dominant  frequencies,  damping 
and  excitation  amplitudes.  Within  the  accuracy  of  the 
measurements,  the  presence  of  the  phot on- Induced  electron 
cloud  had  no  effect  on  the  period  or  damping  tines  of  the 
RESMOD  modes.  The  proximity  of  an  electrical  pulsar  and 
Its  200-0hm  terminating  resistor  markedly  decreaced  the 
damping  time  during  electrical  excitation. 

Prony,  R. ,  "Essal  Experimental  et  Analytlque  sur  les  Lois  de  la  Dllatabllite 
de  Fluldes  Elastlques  et  sur  Celles  de  la  Force  Expansive  de  la  Vaperu  de 
I'Alkool,  a  Differences  Temperatures,"  J.  I'Ecole  Polvtech.  (Paris), 

Vol.  1,  No.  2,  pp.  24-76,  1795. 

This  paper  has  been  abstracted  in  the  Section  2.0  of  this 
bibliography.  This  is  Prony's  original  paper. 


Quasi,  A.H.  and  F.R.  Splcznogla,  "Time-Limited  Signal  Representation  and 
Analysis  by  Using  Non-Fourier  Techniques,"  U.S.  Navy  Journal  of  Under¬ 
water  Acoustics,  Vol.  18,  No.  3,  pp.  397-402,  July  1968. 

Tlme-llffllted  signals  are  represented  and  analyzed  by  using 
non-Fourier  techniques.  Results  show  that  non-Fourier 
techniques  can  be  superior  to  the  conventional  Fourier 
approach  for  the  frequency  analysis  of  time  limited  signals. 


Salz«r,  H.E. ,  "Some  Extensions  o£  Prony  Approximation",  Zeltachrlft  fur 
Angewa^te  Mathematlk,  und  Mechanlk.  57,  /M,  pp.  269-271,  1977. 


Prony *8  method  of  approximating  a  function  f(8)  by  a  sum 

n  a^x 

of  exponential  terma  A^  e  ^  ,  where  both  a^  and  A^ 

are  to  be  determined,  la  extended  to  oscillatory  Interpola¬ 
tion  of  any  order,  and  also  to  the  direct  global  approxima¬ 
tion  to  the  solution  of  the  linear  differential  equation 

£^.0  ^r^*^  -  f(x)  ^0. 

Sarkar,  T.K. ,  J.  Nebat,  and  D.D.  Weiner,  "Suboptlmal  System  Approximation/ 
Identification  with  Known  Error",  Air  Force  Weapons  Laboratory,  AFWL-TR- 
77-200,  September  1977.  Also  Mathematics  Notes  (1149* 

This  paper  presents  a  noniterative  method  for  approxlmat;Lng 
empirical  signals  over  [O,*]  by  a  linear  combination  of 
exponentials.  The  technique  results  In  a  suboptlmal 
approximation.  Notably,  the  dependence  of  the  suboptlmal 
exponents  sj^  on  the  Integral  square  error  c  Is  such  that 

11m  (c-*-)  sj^  <■  s^,  the  optimal  exponents.  The  method  may 

also  be  used  for  system  Identification.  It  Is  especially 
useful  when  the  system  Is  modelled  by  a  black  box  and  one 
has  access  only  to  the  Input  and  output  terminals  of  the 
system.  A  technique  la  demonstrated  to  find  the  multiple 
poles  of  a  system  along  with  the  residues  at  the  poles  when 
the  output  of  the  system  to  a  known  Input  is  given.  Among 
the  advantages  of  the  method  are  its  natural  insensitivity 
to  noise  In  the  data  and  the  explicit  determination  of  the 
signal  order.  Representative  computations  are  made  of  the 
poles  from  the  transient  response  of  a  conducting  pipe 
tested  at  the  ATHAMAS-1  EMF  simulator. 

Sarkar,  T.K. ,  V.K.  Jain,  J.  Nebat,  and  D.D.  Weiner,  "A  Comparison  of  the 
Pencll-of-Functlon  Method  with  Prony's  Method,  Weiner  Filters  and  Other 
Identification  Techniques",  Rochester  Institute  of  Technology,  AFWL 
Mathematics  Notes  /i(54,  December  1977. 

The  problem  of  Interest  Is  to  Identify  the  transfer 
function  of  a  system  by  Its  poles  and  residues  when  the 
noise  contaminated  Input  and  output  are  specd.fied.  The 
first  aim  of  this  paper  Is  to  Illustrate  that  several 
different  formulations  for  characterizing  the  Impulse 
response  of  a  system  yield  the  same  set  of  poles  as  Is 
obtained  In  the  case  of  a  discrete  Weiner  filter.  The 
second  goal  Is  to  show  how  different  formulations 
regularize  the  Ill-posed  system  Identification  problem. 


It  is  dsmonstratad  that  the  Walnar  flltar  la  not  always 
reallzabla  as  a  causal  rational  function.  Whan  tha  ordar 
of  the  flltar  Is  specified  a  priori,  the  resulting  filter 
may  no  longer  be  optimum.  Finally,  representative 
computations  are  made  of  tha  poles  from  the  transient 
response  of  a  conducting  pipe  tested  at  the  ATHAMAS-1  EMP 
simulator  to  demonstrate  the  stability,  reliability, 
consistency,  and  accuracy  of  the  results  obtained  by  the 
pencilz-of-f unction  method. 

Sarkar,  T.K. ,  D.D.  Weiner,  and  V.K.  Jain,  "Some  Mathematical  Considerations 
In  Dealing  with  the  Inverse  Problem",  IEEE  Transactions  on  Antennas  and 
Propagation.  Vol.  AP-29,  No.  2,  pp.  373-379,  March  1981. 


Many  problems  of  mathematical  physica  can  be  formulated  In 
terms  of  tha  operator  equation  Ax  ■  y,  where  A  Is  an 
Integro-dlfferential  operator.  Given  A  and  x,  the  solution 
for  y  is  usually  straightforward.  However,  the  inverse 
problem  which  consists  of  tha  solution  for  x  when  given  A 
and  y  is  much  more  difficult.  The  following  questions 
relative  to  the  Inverse  problem  are  explored.  1)  Does 
specification  of  the  operator  A  determine  the  set  {y}  for 
which  a  solution  x  Is  possible?  2)  Does  the  Inverse  problem 
always  have  a  unique  solution?  3)  Do  small  perturbations 
of  tha  forcing  function  y  always  result  In  small  perturba¬ 
tions  of  the  solution?  4)  What  are  some  of  the  considera¬ 
tions  that  enter  Into  the  choice  of  a  solution  technique 
for  a  specific  problem?  The  concept  of  an  ill-posed 
problem  versus  that  of  a  wall-posed  problem  Is  discussed. 
Specifically,  the  manner  by  which  an  Ill-posed  problem  may 
be  regularized  to  a  well-posed  problem  is  presented.  The 
concepts  are  illustrated  by  several  examples. 

Scharf,  L.L.  and  A. A.  Beex,  "Covariance  Sequence  Approximation  for 

Parametric  Spectrum  Modelling",  Technical  Report  August  80-1,  Colorado 

State  University,  Ft.  Collins,  Colorado  80523. 


Parametric  methods  of  spectrum  analysis  are  founded  on 
flnlte-dlmenslonal  models  for  convarlance  sequences. 
Rational  spectrum  approxlmants  for  continuous  spectra  are 
based  on  autoregressive  (AR) ,  moving  average  (MA)  or  auto¬ 
regressive  moving  average  (ARMA)  models  for  covarlauce 
sequences.  Line  spectrum  approxlmants  to  discrete  spectra 
are  based  on  cosinusoidal  models  for  covariance  sequences. 
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In  this  paper  the  authors  make  the  point  that  a  wide  CJ 

variety  of  spectrum  types  admit  to  modal  analysis  wherein 

the  modes  are  characterized  by  amplitudes,  frequencies, 

and  damping  factors.  The  associated  modal  decomposition 

Is  appropriate  for  both  continuous  and  discrete  components 
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of  cha  spactrutn.  ...a  domain  of  attraction  for  tha 
dacompoaltlon  Includas  ABMA  aaquancaa*  harmonically-  or 
nonharmonically-rclated  alnuaolda,  danspad  slnuaolda,  whlta 
noise,  and  lluaar  combinations  of  these. 

fha  parametric  spectrum  analysis  problem  now  becomes  one  of 
identifying  mode  parameters.  This  is  achieved  by  solving 
two  modified  least-squares  problems,  tlumarical  results  are 
presented  to  illustrate  the  identification  of  mode 
parameters  and  corresponding  spectra  from  finite  records  of 
perfect  and  estimatad  covariance  sequences.  The  results 
for  sinusoids  and  sinusoids  in  white  noise  are  interpreted 
in  terms  of  inphase  and  quadrature  effects  attributable  to 
the  finite  record  length. 

Scharf,  L.L.,  A. A.  Beex  and  T.  von  Reyn,  "Modal  Decomposition  of 
Covariance  Sequences  for  Paranwtrlc  Spectrum  Analysis",  contributed  peper, 
1981  International  Conference  on  Acoustics,  Speech,  and  Signal  Processing, 
Atlanta,  March  30,  1981. 

In  this  paper  the  authroa  make  the  point  that  a  wide 
variety  of  spectrum  types  admit  to  modal  analysis  wherein 
the  modes  are  characterized  by  amplitudes,  frequencies, 
and  damping  factors.  The  associated  modal  decomposition  is 
appropriate  for  both  continuous  and  discrete  components  of 
the  spectrum.  The  domain  of  attraction  for  the  decomposi¬ 
tion  Includes  ARMA  sequences,  harmonically-  or  nonharmonl- 
cally-related  sinusoids,  damped  sinusoids,  white  noise,  end 
linear  combinations  of  these. 

Numerical  results  are  presented  to  illustrate  the  Identifi¬ 
cation  of  mode  parameters  and  corresponding  spectra  from 
finite  records  of  perfect  and  estimated  covariance  sequences. 

The  results  for  sinusoids  and  sinusoids  in  white  noise  are 
interpreted  in  terms  of  inphase  and  quadrature  effects 
attributable  to  the  finite  record  length. 

Schaubert,  D.H. ,  "Application  of  Prony's  Method  to  Time-Domain 
Ref lactometer  Data  and  Equivalent  Circuit  Synthesis",  HDL-TR-1857,  Harry 
Diamond  Laboratories,  June  1978.  Also  IEEE  Transactions  on  Antennas  and 
Propagation.  Vol.  AP-27,  No.  2,  pp.  180-184,  March  1979. 

The  singularity  expansion  of  time-domain  reflectometer  data 
has  been  computed  and  used  to  calculate  lumped-element 
equivalent  circuits  for  the  impedance  of  some  typical 
antennas.  Prony's  algorithm  was  used  to  obtain  the  poles 
of  the  artennas'  terminal  voltage  veveform  due  to  a  step¬ 
like  excitation.  The  residues,  however,  were  calculated 
subject  to  the  uniformrsrror  uorm  Instsed  of  the  leest- 
squarsd-error  norm.  Although  it  requires  much  longer  to 
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calculaca,  th«  unlform-arror  norm  approximation  la  of tan 
orders  of  magnitude  better  than  the  laast-squared-error 
approximation.  The  phyalcal  realizability  of  Impedance 
functions  obtained  from  experimental  data  has  been 
Investigated.  AlsOi  the  procedures  necessary  to  synthesize 
the  lumped-element  network  have  been  evaluated.  A  simple 
partial-fraction  expansion  of  a  realizable  Impedance 
function  representing  an  antenna  was  found  to  yield 
nonreallzabla  element  values.  A  general  network  synthesis 
procedure,  such  as  Brune's  method,  Is  required. 

Schultz,  B.R. ,  "Estimation  of  Pulse  Transfer  Function  Paramaters  by  Qua. 

linearization",  IEEE  Transactions  on  Automatic  Control,  pp.  424-A26, 

August  1963. 

An  Iterative  algorithm  for  the  Identification  of  single¬ 
input  single-output  linear  stationary  discrete  system  is 
developed  using  the  method  of  quaalllnearlzatlon.  The 
resulting  procedure  Is  similar  to  mode  1  of  the  method  of 
Stelglltz  and  McBride  but  has  the  advantage  of  the 
quadratic  convergence  property  of  quasillnearlzatlon.  It 
la  shown  that  this  algoirlthm  becomes  mode  1  If  the 
measured  plant  output  Is  used  in  the  calculations  In  place 
of  the  model  output.  Consequently,  the  two  methods  are 
extremely  compatible  and  It  is  a  simple  matter  to  combine 
them  In  a  single  program,  which  generates  Its  own  Initial 
estimate,  has  the  wide  range  of  convergence  of  mode  1,  and 
possesses  the  quadratic  convergence  property  of  quaal- 
linearlzatlon  for  final  convergence  to  a  solution.  The 
method  also  permits  the  estimation  of  plant  Initial 
conditions  In  those  cases  where  they  must  be  considered. 

Results  of  a  few  numerical  applications  are  discussed. 

Scrlvner,  G.J.,  "Prony's  Analysis  in  the  Presence  of  Noise",  Computer 

Sciences  Corporation,  Report  II/3-C-23S,  October  1977.  Also  AFWL 

Mathematical  Note  i!l67. 

This  report  studies  a  methodology,  based  on  Prony's 
algorithm,  for  extracting  complex  frequencies  and 
associated  residues  directly  from  transient  response  data. 

The  methodology  is  not  new  but  has  been  receiving 
Increasing  attention  due  to  its  potential  usefulness  In  a 
variety  of  different  applications  In  quite  diverse  scien¬ 
tific  disciplines.  Although  this  particular  study  was  to 
evaluate  the  potential  usefulness  of  Prony's  method  for 
analysis  of  System  Generated  Electromagnetic  Pulse  (SGEMP) 
experimental  data,  the  results  obtained  are  of  a  much 
broader  nature.  The  report  presents  a  classical  derivation 
of  the  matrix  aquation  which  result  if  one  uses  Prony's 
algorithm  in  conjunction  with  a  least-squares  criterion. 


Tha  resulting  matrix  Is  atudlad  with  rsspact  to  tho  nature 
of  Its  alganvalua  structure  and  an  error  analysis  Is 
developed  which  highlights  the  Importance  of  this  structure. 

The  influence  of  noise  on  the  eigenvalue  structure  and  Its 
impact  on  the  determination  of  the  number  of  poles  in  the 
data  are  examined.  Optimisation  concepts  relative  to 
advantageous  modification  of  the  eigenvalue  structure  of 
the  Prony  matrix  ''re  qualitatively  discussed.  Difficulties 
associated  with  the  two  state  application  of  the  least- 
squares  method  In  conventional  Prony  analysis  are  cited  and 
an  iterative  method  of  removing  this  shortcoming  described. 

Numerical  calculations  illustrating  the  majority  of  the 
above  concepts  are  presented  and  discussed. 

Siegmau,  E.  and  H.Y.  Miller^  "Unstable  Optical  Resonator  Loss  Calculations 
I  Using  the  Prony  Method",  Applied  Optics.  Vol.  9,  No.  12,  pp.  2729-2736, 

December  1970. 

The  eigenvalues  for  all  the  significant  low-order  resonant 
modes  of  an  unstable  optical  resonator  with  clrcilar 
mirrors  are  computed  using  an  eigenvalue  laethod  called  the 
)  Prony  method.  A  general  equi valance  relation  is  also 

given,  by  means  of  which  one  can  obtain  tha  design 
parameters  for  a  single-ended  unstable  resonator  of  the 
type  usually  employed  In  practical  lasers,  from  the 
calculated  or  tabulated  values  for  an  equivalent  symmetric 
or  double-ended  unstable  resonator. 

I 

Simon,  W. ,  "A  Method  of  Exponential  Separation  Applicable  to  Small 
Computers",  Phvs.  Med.  Biol..  Vol.  15,  No.  2,  pp.  355-360,  1970. 

A  method  has  been  developed  for  separation  of  the  components 
of  multlcompartmsnt  exponentials  which  operates  on  a  small 
I  computer.  It  seems  free  of  many  of  the  difficulties  of 

ocher  methods  and  provides  a  reasonable  test  for  the  number 
of  components. 

Smith,  W.D.  and  D.L.  Lager,  "Parametric  Characterization  of  Random 
Processes  Using  Prony* s  Method",  Lawrence  Livermore  Laboratory,  UCRL- 
>  52673,  25  May  1979. 

The  compact  parametric  characterization  of  a  random 
process  Is  often  valuable  in  problems  involving  signal 
classification  and  system  identification.  Hie  method  of 
Prony  suggest  two  different  approaches  to  obtaining  such  a 
*  characterization  from  a  finite  record  of  sampled  data. 

The  first  approach  is  to  extract  the  parameters  of  the 
estimated  autocorrelation  waveform.  The  second  approach  is 
to  estimate  Che  parameters  of  a  white  noise  driven  linear 
system  that  generates  an  output  with  Che  same  character¬ 
istics  as  the  given  random  process. 


Thast  two  approaches  are  compared  with  ocher  common  signal 
processing  mechodo.  Also,  Che  performance  of  each  approach 
is  evaluau-'d  cheoracically  and  wich  compucer^slmulaced  and 
experi^uncally  recorded  data.  Each  approach  is  found  co 
rapidly  characterize  both  a  ationary  and  nonstationary  random 
proceases.  The  first  approach  is  superior  for  the  accurate 
paramatric  charactevlzatiou  of  a  stationary  process.  Alter¬ 
nately,  the  more  precise  second  approach  is  best  for  monitoring 
changes  vith  time  in  random  process  parameters. 

Soussou,  J.E.,  F.  Moavanzadah,  and  M.H.  Gradowczyk,  "Application  of  Prony 
Series  to  Linear  Vsicoelastlclty",  Transactions  of  the  Society  of 
Rheology.  Vol.  14,  No.  4,  pp.  573-584,  1970. 

A  new  method  for  description  of  viscoelastic  functions  with 
Prony  series  of  exponentials  is  presented  and  compared  with 
existing  methods.  This  method  pieces  constraints  on  the 
coefficients  of  the  series  which  guarantees  smoothness  of 
thu  functions  and  a  discrete  spectral  representation.  An 
optimization  technique  is  used  to  obtain  the  solution.  The 
resulting  series  permit  the  use  of  Whittaker's  method  of 
solution  of  integral  equations  which  is  specialized  to  the 
use  of  viscoelastic  analysis,  and  an  exact  relauionship  is 
obtained  between  creep  and  relaxation  functions. 

Spitznogle,  F.R.  and  A.H.  Quazi,  "Representation  and  Analysis  of  Tiue- 
Liffldted  Signals  Using  a  Complex  Exponential  Algorithm",  Journal  of  the 
Acoustical  Society  of  America.  Vol.  47,  No.  5,  pp.  1150-1155,  May  1970. 

A  complex  exponential  algorithm  developed  for  the 
representation  and  analysis  of  time-limited  signals  is 
defined,  and  its  evaluation  with  respect  to  conventional 
discrete  Fourier  techniques  is  discussed.  It  is  shown  that 
for  a  given  length  of  a  signal  containing  discrete  frequency 
information,  sampled  at  least  to  the  Nyqulst  criterion,  the 
complex  exponeatlal  algorithm  can  often  provide  increased 
frequency  resolution  over  standard  Fourier  techniques.  It 
is  also  shown  that  the  complex  exponential  algorithm 
provides  an  improved  mechanism  over  Fourier  techniques  for 
interpolation  between  points  in  a  sampled  signal  containing 
discrete  frequency  components  in  the  presence  of  board-band 
noise.  The  effects  of  noise  in  the  complex  exponential 
technique  and  the  computational  difficulties  associated  with 
the  present  complex  exponential  algorithm  are  also  discussed. 


Texas  Inatrumente  Incorporated,  "Annual  Report  -  Representation  and 
Analysis  of  Sonar  Signals  Volume  I:  Improvements  In  the  Complex  Ex¬ 
ponential  Signal  Analysis  Computational  Algorithm,"  Texas  Instruments, 
Inc.,  Dallas,  Tr  *.as,  31  March  1970. 

This  volume  of  the  1969-70  Annual  Report  on  Contract 
N00014-69-C-031S  describes  the  results  of  a  study  to  re¬ 
duce  the  computer  memory  and  computational  time  required 
to  obtain  a  complex  exponential  signal  representation. 

As  described  In  the  report,  substantial  reductions  In  com¬ 
putational  complexity  have  been  achieved  through  improvements 
In  the  Prony  method  for  solving  the  nonlinear  system  of 
equations  associated  with  a  complex  exponential  signal 
analysis.  Projections  of  the  real-time  capability  of  a 
hard-wired  complex  exponential  processor  using  state-of- 
the-art  hardware  are  given.  In  addition,  two  Iterative 
complex  exponential  algorithms  based  on  the  Mewton-Raphson 
and  gradient  techniques  for  solving  a  nonlinear  system  of 
equations  are  described. 

Trlvett,  David  H.,  "A  Modified  Prony  Method  Approach  to  Echo-Reduction 
Measurements  of  Time-Limited  Transient  Signals,"  MRL  Memorandum  Report 
4172,  June  16,  1980,  ADA  087211. 

A  modified  Prony  method  Is  presented  for  measuring  the 
steady-state  echo  reduction  of  acoustic  panels.  The 
method  extrapolates  the  stead-state  amplitudes  from  the 
transient  portion  of  the  signal  allowing  time-limited 
measurements.  The  method  is  applied  to  measurements  of 
square  panels  76  cm  on  an  edge  and  0.95  cm- thick  steel 
and  aluminum  in  the  frequency  range  of  3  -  10  kHz.  The 
signals  were  time-limited  to  200ys  (0.6  -  2»0X)  by  the 
arrival  of  the  diffracted  signal  from  the  panel  edges. 

Results  are  compared  with  theoretical  values  and  a  pro¬ 
gram  listing  Is  included. 


Van  Blarlcum,  M.L.  and  R.  Mlttra,  "A  Technlr|ue  for  Extracting  the  Poles 
and  Roaldues  of  a  System  Directly  from  Its  Transient  Response,"  IEEE 
Transactions  on  Antennas  and  Propagation.  Vol.  AP-23,  No.  6,  pp.  777- 
781,  November  1975. 


A  novel  approach  for  systematically  deriving  the  complex 
poles  and  residues  of  a  structure  from  a  set  of  time-domain 
data  is  presented.  The  method  Is  based  on  Prony's  algorltluo 
which  Involves  the  Inversi.a  of  two  matrices  and  a  solution 
of  the  zeros  of  an  nth  degree  polynomial,  n  being  the  number 
of  desired  poles.  Two  numerical  examples  are  presented  and 
several  numerical  advantages  of  this  method  are  discussed. 

Van  Blarlcum,  M.L.,  "Techniques  for  Extracting  the  Complex  Resonances  of 
a  System  Directly  from  Its  Transient  Response,"  Ph.D.  Dissertation,  De¬ 
partment  of  Electrical  Engineering,  University  of  Illinois,  December  1975; 
also  AF^^  EMP-Interactlon  Note  3ri. 


The  singularity  expansion  method  Is  a  techniques  used  for 
writing  the  transient  response  of  a  structure  as  a  sum  of 
exponentially  damped  sinusoids.  In  order  to  apply  this 
method.  It  Is  first  necessary  to  obtain  the  singularities 
of  the  system  being  studied.  The  conventional  approach  for 
determining  the  singularities  of  a  system  Is  based  on  an 
iterative  search  procedure  that  seeks  the  zeros  of  the 
system  determinant  In  the  complex  frequency  plane.  The 
alternative  approach  of  extracting  the  system  singularities 
directly  t  rom  the  transient  response  function  Is  discussed. 
The  method  developed  Is  based  on  Prony's  algorithm.  Some 
basic  problems  which  are  associated  with  the  use  of  Prony's 
algorithm  are  discussed  and  solutions  are  obtained. 


It  is  demonstrated  that  Prony's  method  Is  appllcablo  to 
systems  with  multiple  as  well  as  simple  pole  singularities. 
Two  techniques  are  presented  for  systematically,  determining 
the  number  of  poles  contained  In  a  transient  response. 


Van  Blarlcum,  M.L.  and  R.  Mlttra,  "FSPL:  A  Computer  Program  for 
Extracting  the  Poles  and  Residues  from  a  Transient  Signal  Using  Prony's 
Method",  Electromagnetics  Laboratory  Scientific  Report  No.  7602, 
University  of  Illinois,  Urbana,  Illinois,  March  1976. 

This  report  describes  the  Fortran  program  PSPL,  which  Is 
used  to  extract  the  poles  and  residues  of  a  system  from  a 
set  of  discrete  transient  data  using  a  least-squares 
version  of  Prony's  method.  The  code  has  been  written  so 
chat  multiple  as  well  as  simple  poles  can  be  obtained  If 
they  are  present  In  the  data.  The  listing  of  Che  program 
given  In  this  report  Is  for  use  on  the  CDC-7600  computer 
but  can  be  easily  adapted  to  other  machines. 
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Van  Blarlcua,  M.L.  and  R.  Mlttra,  "FEIGEN:  A  Computer  Program  for  the 
Determination  of  the  Number  of  Poles  In  a  Transient  Signal", 

Electromagnetics  Laboratory  Report  No.  76-5,  University  of  Illinois, 

Urbana,  Illinois,  June  1976. 

This  la  a  users  program  for  the  PEIGEN  program  which  Is 
used  to  determine  the  number  of  poles,  N  In  a  set  of 
transient  data  assuming  that  the  data  are  made  up  of  a 
finite  sum  of  exponentials.  PEIGEN  Is  written  as  a 
companion  to  the  program  PSPL,  which  Is  used  to  extract  the 
N  poles  from  the  transient  data  once  the  value  of  N  has 
been  determined.  The  method  uses  the  eigenvalue  method. 

Van  Blarlcum,  M.L. ,  "The  Use  of  Prony’s  Method  as  an  Analysis  Tool  Applied 
to  Electrical  Excitation  of  FltSatCom",  Mission  Research  Corporation  Report, 
MRC-R-3C\,  January  1977. 

This  report  discusses  the  use  of  the  numerical  algorithm 
known  as  Prony's  method  as  a  tool  for  the  analysis  of 
electrical  transient  experimental  data.  The  data  to  which 
these  techniques  are  applied  as  examples  are  results  of 
electrical  excitation  tests  of  a  structural  mockup  of  the 
FltSatCom  satellite. 

Van  Blarlcum,  M.L.  and  R.  Mlttra,  "Problems  and  Solutions  Associated  with 
Prony's  Method  for  Processing  Transient  Data",  IEEE  Transactions  on 
Antennas  and  Propagation.  Vol.  AP-26,  No,  1,  pp,  174-182,  January  1978, 

Three  difficulties  associated  with  Prony's  method  are 
studied.  These  are  the  extension  of  the  method  to  allow 
multiple  poles,  the  development  of  techniques  for 
determining  the  number  of  poles  contained  In  the  transient 
data,  and  the  effects  of  noise  In  the  data  on  the  numerical 
procedures.  Solutions  to  these  difficulties  are  studied, 
and  numerical  samples  are  presented. 
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Weaver,  J.J.,  D.  Ng,  and  D.  L.  Lager,  "Analytical  Qualification  of  System 
Identification  (Mod^  Analysis)  Codes  for  Use  In  the  Dynamic  Testing  of 
Nuclear  Power  Plant  Structure",  Lawrence  Livermore  Laboratory,  UCID- 18144, 
2  January  1980. 

This  report  describes  the  analytical  evaluation  of  two 
particular  system  Identification  codes  used  at  Lawrence 
Livermore  Laboratory.  Both  codes  are  elgenparameter  Iden¬ 
tification  codes,  however,  one  uses  a  time  domain  approach 
while  the  other  a  frequency  domain  approach.  The  evaluation 
was  accomplished  by  analytically  generating  several  time 
history  signals  In  which  the  true  modal  parameters  were  known. 

These  time  histories  ranged  from  widely  spaced  modes  with 
spacing  factors  of  6  percent.  These  signals  were  then 
polluted  with  vai.'loua  levels  of  simulated  measurement  noise 
and  the  ability  of  our  computer  codes  to  extract  the  para¬ 
meters  from  this  noisy  data  was  evaluated. 

Weiss,  L.  and  R.M.  McDonough,  "Prony's  Method,  Z-Transforms ,  and  Pade 
Approximation",  SIAM  Review.  Vol.  5,  No.  2,  pp.  145-149,  April  1963. 

The  purpose  in  this  note  Is  to  point  out  that  the  original 
form  of  Prony's  method  may  also  be  regarded  as  a  Pade 
approximation  In  the  Z-transform  domain. 

Weygandt,  C.M.  and  N.N.  Purl,  "Identification  of  a  Linear  System  From 
Discrete  Values  of  Input-Output  Data",  1966  Proceeding  4th  Annual 
Allerton  Conference  on  Circuit  and  System  Theory,  po.  903-912,  1966. 

In  this  paper  the  authors  describe  a  method  for 
determining  the  transfer  function  of  a  physical  system 
having  period  Input  and  output,  both  of  which  can  be 
measured.  The  method  depends  upon  taking  2n  equidistant 
samples  of  the  Input  and  the  output  within  one  period,  n 
being  the  number  of  poles  the  system  transfer  function  Is 
assumed  to  have.  These  samples  are  used  to  compute  (n+l) 
determinants  and  thereby  the  coefficients  and  the  roots  of 
an  nth  order  polynomial.  The  roots  of  this  nth  order 
polynomial  are  further  processed  to  obtain  the  poles  of  the 
transfer  function.  Once  the  location  of  poles  Is 
determined,  a  very  simple  computer  routine  determines  the 
residues  of  these  poles  and  thereby  the  zeros  and  the  gain 
constant  are  determined.  An  essential  feature  of  the  scheme 
Is  that  it  determines  the  poles  Independent  of  the  zeros. 

In  practical  examples  the  number  of  samples  taken  IS  much 
larger  than  2n,  thereby  generating  a  large  number  of  nth 
order  polynomials.  The  average  of  the  roots  of  these 
polynomials  are  used  to  locate  the  poles. 
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Vftilte,  O.R.,  "Inversion  of  the  Llmb-Darkenlng  Equation  Using  the  Frony 
Algorithm",  The  Astophyslcal  Journal.  Vol.  152,  pp.  217-223,  April  1968. 

King's  Interesting  application  of  the  Frony  algorithm  to 
the  Inversion  of  the  llmb-darkenlng  equation  Is  modified 
and  generalized  so  that  It  can  be  applied  rigorously  to 
llmb-darksnlng  curves,  spectral  line  profiles,  and 
multlplet  line  Intensities.  Since  King's  formulation 
requires  physically  Inaccessible  observations,  a  change  of 
variable  Is  Introduced  In  the  llmb-darkenlng  Integral  to 
avoid  this  difficulty.  Calculations  on  noisy  data 
confirm  the  self-llmltlng  property  of  this  Inversion 
method.  An  approximate  method  for  constructing  the  smooth 
S(t)  distribution  from  the  slab  solution  Is  presented. 

Wlscombe,  W.J.  and  J.  W.  Evans,  "Exponential-Sum  Fitting  of  Radiative 
Transmission  Functions",  Journal  of  Computational  Fhvslcs.  24,  pp.  416- 
444,  1977. 

Fitting  transmission  functions  with  exponential  sums  Is  the 
basis  for  a  widely  used  approximation  for  calculating 
spectrally  Integrated  radiative  fluxes  In  planetary 
atmospheres,  expeclally  when  both  line  absorption  and 
scattering  are  Important.  The  error  In  this  method 
depends  crucially  on  the  accuracy  of  the  fits,  but 
unfortunately  exponential-sum  fitting  Is  a  classical  Ill- 
conditioned  problem  of  numerical  analysis.  Frevlous 
techniques  devised  for  exponential- sum  fitting  are  often 
unsatisfactory  In  this  application.  This  part  presents  a 
new  method  which  sidesteps  the  111-condltloning,  guarantees 
convergence  to  the  unique  best  least-squares  fit,  gives 
positive  coefficients,  and  produces  fits  orders  of 
magnitude  more  accurate  than  any  which  have  so  far  been 
published.  The  method  Is  demonstrated  to  be  capable  of 
recovering  an  exponential  sum,  given  data  sampled  from  that 
sum  and  rounded  to  as  few  os  two  decimal  places.  Sample 
fits  ere  given  for  the  Goody  and  Malkmus  random  band 
models  and  the  Yamamoto  H2O  solar  absorption  data  In  order 
to  Illustrate  the  high  accuracy  of  the  method.  The  effect 
on  the  fits  or  errors  In  the  transmission  data  Is  examined 
in  some  depth. 

Wolf,  J.K. ,  "Decoding  of  Bose-Chaudhurl-Hocquenghem  Codes  and  Frony 's 
Method  for  Curve  Fitting",  IEEE  Transactions  on  Information  Theory,  p.  608, 
October  1967. 

The  literature  Is  filled  with  examples  of  a  common  set  of 
equations  which  arise  in  two  or  more  diverse  applications. 

The  purpose  of  this  correspondence  Is  to  point  out  that 
such  a  situation  has  occurred  In  the  two  fields  of 


1)  algebraic  coding  theory,  and  2)  curve  fitting.  Two 
aspects  are  of  partlcxilar  Interest:  the  large  time  span 
which  separates  the  fundamental  works  In  these  two  areas, 
and  that  the  methods  of  solution  are  Identical. 

Wynn,  P. ,  "A  Note  on  the  Fitting  of  Certain  Types  of  Experimental  Data", 
Statistlca  Neerlandlca.  Vol,  16,  No.  2,  pp.  143-150,  1962. 

This  paper  Is  concerned  with  the  problem  of  determining 
those  values  of  the  parameters  Z,  Ai,  ai(i-l,2, . . . ,h) 
which  provide  a  best  fit  of  the  form 

k  -a .  t 

f*(t)  -  Z  +  E  A.e  ^ 
i-1 

to  the  obser'ed  values  fg  ■  f (tg) (s*l,2, . . . ,n)  in  the  sense 

that  the  s-jun  E  {t*(r.g)  -  f(tg)}  is  a  minimum. 
i-1 

Two  iterative  numerical  techniques  for  the  solution  of  this 
problem  are  discussed.  The  first,  a  variable  gradient  method, 
converges  rapidly  but  involves  a  relatively  large  amount  of 
computation;  the  second  involves  less  computation  but  con¬ 
verges  less  rapidly.  It  is  shown  that  by  application  of  a 
suitable  acceleration  technique  to  the  second  method,  its 
con  'ergence  is  made  as  rapid  as  that  of  the  first. 
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Yengsc,  W.C.,  "Approximation  to  a  Specified  Time  Response,"  IRE  Trans¬ 
actions  on  Circuit  Theory,  pp.  152-163,  June  1962. 

This  paper  presents  a  procedure  by  which  specified  data  or  •\ 
function  of  time  h^Ct)  can  be  approximated  by  trigonometric 
and/or  exponential  functions  of  time  h(t)  for  which  the 
Laplace  transformations  H(s)  are  known  and  can  be  expressed 
In  rational  fraction  form.  The  procedure  Is  based  on  fitting 
h((t)  by  an  mth-order  difference  equation  whose  coefficients 
are  determined  by  a  least-squares  technique.  These  co¬ 
efficients  are  used  directly  to  determine  the  poles  of  H(s). 

The  zeros  of  H(s)  are  established  by  using  the  prescribed 
data  or  function  h*(t)  and  the  Initial  value  theorem.  The 
approximate  function  of  time  is  obtained  by  taking  the  In¬ 
verse  Laplace  transformation  of  H(s) .  By  this  procedure 
not  only  Is  an  approximation  obtained  for  h*(t)  In  the  time 
domain,  but  Its  transform  la  also  found  In  rational  fraction 
form  suitable  for  realization  as  a  driving  point  or  transfer 
function.  Furthermore,  the  least-squares  technique  used  In 
determining  moat  or  all  of  the  unknown  parameterN  In  this 
procedure  tends  to  minimize  the  effect  of  random  errors  or 
noise  present  in  the  specified  data. 

Youla,  D.C.,  "The  Identification  of  Linear  Dynamical  Systems  from  Time- 
Domain  Measurements,"  Department  of  Electrical  Engineering  -  Farmlngdale, 
New  York  11735,  pp.  1-26,  Appendix  A1-A2,  1  September  1978. 

The  problem  of  Identifying  the  impulse  response  of  a  lumped, 
linear  time- Invariant  system  from  a  discrete  set  of  Its 
samples,  Is  essentially  flnlte-dlmenslonal  In  character. 

Frony's  method  is  distinguished  by  t^a  fact  that  It  re¬ 
cognizes  this  finite-dimensionality  from  the  outset  and 
thereby  succeeds  in  avoiding  otuch  of  the  numerical  In¬ 
stability  exhibited  by  algorithms  that  are  designed  to 
work  In  an  infinite-dimensional  setting.  This  paper 
attempts  to  generate  some  new  Insight  into  both  the  theo¬ 
retical  and  practical  implications  of  Frony's  ideas. 

Young,  T.Y.,  "Representation  and  Analysis  of  Signals,  Part  XVI,  Repre¬ 
sentation  and  Detection  of  Multiple-Epoch  Signals,"  The  Johns  Hopkins 
University,  Baltimore,  Maryland,  AS-A05028,  May  1963. 

A  multiple-epoch  signal  consists  of  several  signals  which 
for  one  reason  or  another,  overlap  one  another.  The  purpose 
of  this  report  Is  to  desig^n  a  procedure  for  the  detection 
of  the  Individual  epochs  of  the  overlapping  signals  and  to 
represent  them  properly.  The  multiple-epoch  signal  may  be 
corrupted  with  random  Gaussian  noise  of  zero  mean.  The 
individual  signals  are  assumed  to  be  representable  by  a 
set  of  exponential  functions  with  acceptable  error,  and  the 
exponents  of  this  set  of  exponentials  are  assumed  known. 

Any  two  adjacent  epochs  are  assumed  to  be  separated  at 
least  Tg  seconds  apart. 

For  the  signal  uncorrupted  with  noise  a  criterion  based  on 
the  error  energy  is  described  which  is  of  theoretical  interest. 


This  crltsrlon  Is  useful  for  the  detection  of  first  epoch 
only.  To  detect  the  second  epoch.  Che  first  slgnel  Is  sub¬ 
tracted  from  the  original  signal  by  means  of  the  "comple¬ 
mentary*  operator  concept.  For  a  noisy  signal,  a  likeli¬ 
hood  ratio  criterion  Is  proposed.  The  preliminary  experi¬ 
mental  results  carried  out  on  a  digital  computer  Justify  the 
theoretical  study. 

Young,  T.Y.  and  W.H.  Huggins,  "Representation  and  Analysis  of  Signals 
Part  XIII:  Discrete  Orthonormal  Exponentials,"  The  Johns  Hopkins  Univer¬ 
sity,  Department  of  Electrical  Engineering,  Baltimore,  Maryland,  AS- 
405028,  May  1963. 

In  dealing  with  the  problem  of  signal  analysis,  one  set  of 
very  useful  component  functions  Is  the  orthonormal  expon¬ 
ential  functions.  Previous  studies  have  concentrated 
mostly,  on  continuous  exponentials  with  real  exponents. 

With  recent  progress  In  high-speed  digital  computer.  It 
Is  desirable  to  perform  the  signal  analysis  on  digital  com¬ 
puter.  This  requires  the  transformation  of  orthonormal 
exponentials  from  the  continuous  form  Into  the  sampled  data 
form,  that  la,  the  z-tranaformatlon. 

This  paper  deals  with  digitizing  the  orthunormal  expo¬ 
nentials  with  complex  exponents.  The  z-transforms  of 
continuous  orthonormal  exponentials  are  not  themselves 
orthogonal  In  z-domaln.  This  Is  due  to  the  fundamental 
property  of  sampling  that  a  signal  with  a  frequency  spectrum 
higher  than  the  sampling  frequency  Is  not  exactly 
reproducible  from  Its  sampled  data.  For  this  reason,  the 
Ksutz  process  Is  Introduced  in  the  z-domain  to  define  a  set 
of  discrete  orthonormal  exponentials.  These  discrete 
orthonormal  exponentials  approach  the  ordinary  continuous 
orthonormal  exponentials  In  the  limit  as  the  sampling 
Interval  T  approaches  zero. 

Young,  T.Y.  and  W.H.  Huggins,  "On  the  Representation  of  Electrocardio¬ 
grams",  IEEE  Transactions  on  Blo-Medlcal  Electronics.  Vol.  BME-10,  No.  3, 
pp.  86-95,  July  1963. 

In  dealing  with  the  sophisticated  statistical  analysis  of 
medical  signals  such  as  the  electrocardiogram  (ECC) ,  the 
first  problem  one  encounters  Is  how  to  describe  each  ECG 
by  a  few  numbers.  This  Is  the  problem  of  efficient  re¬ 
presentation  of  signals,  l.e.,  to  approximate  the  signal 
with  the  smallest  number  of  basis  signals  while  preserving 
the  accuracy  of  the  approximation.  This  paper  begins  with 
a  dlscuHSlon  of  signal  representation  In  general.  The  con¬ 
cept  of  signal  space  Is  introduced,  which  is  very  helpful 
In  understanding  the  Ideas  of  signal  representation.  This 
portion  of  the  material  Is  of  tutorial  nature. 
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Attention  is  then  directed  to  A  set  of  basis  components 
which  have  found  to  be  very  efficient  for  EC6  representa¬ 
tion.  These  components  are  the  so-called  orthonormal 
exponential  signals.  An  Iterative  prricess  Is  developed 
which  enables  one  to  fln'^  a  net  of  mntched  exponents  for 
the  representation  of  el.  '  .  With  six  pairs  of  such 
exi)ouentlals ,  the  average  »  of  ECG  representation  (QRS 
and  T  waves  only,  leaving  ouc  P  wave)  Is  In  the  vicinity 
of  five  percent.  Experimental  results  will  be  shown. 

Utilng  this  representation,  further  statistical  analysis  may 
be  carried  out  with  ease. 
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